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INVESTIGATION OF RADIATIVE INTERACTIONS 
IN SUPERSONIC INTERNAL FLOW 


By 

S.N.Tiwari 1 and A.M. Thomas 2 
SUMMARY 

Analyses and numerical procedures are presented to investigate the radiative 
interactions of absorbing emitting species in chemically reacting supersonic flow in various 
ducts. The two-dimensional time dependent Navier-Stokes equations in conjunction with 
radiative flux equation are used to investigate supersonic flows undergoing finite rate 
chemical reaction in a hydrogen air system. The specific problem considered is the flow of 
premixed radiating gas between parallel plates. Specific attention was directed toward 
investigating the radiative contribution of H 2 0, OH, and NO under realistic physical and 
flow conditions. Results are presented for the radiative flux obtained for different gases and 
for various combination of these gases. The problem of chemically reacting and radiating 
flows was solved for the flow of premixed hydrogen-air through a 10° compression ramp. 
Results demonstrate that the radiative interaction increases with an increase in pressure, 
temperature, amount of participating species, plate spacing, and Mach number. Most of the 
energy, however, is transferred by convection in the flow direction. In general the results 
indicate that radiation can have a significant effect on the entire flow field. 


1 Eminent Professor, Department of Mechanical Engineering and Mechanics, Old Dominion 
University, Norfolk VA 23529. 

2 Graduate Research Assistant, Department of Mechanical Engineering and Mechanics, Old 
Dominion University, Norfolk VA 23529. 
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Chapter 1 


INTRODUCTION 

In the past several years, there has been a great deal of research toward development 
of a hypersonic transatmospheric vehicle. At NASA Langley Research Center, the 
hydrogen-fueled Supersonic combustion ramjet (Scramjet) engine has been a strong 
candidate for propelling such a vehicle. Both experimental and numerical techniques 
are being employed for a better understanding of the complex flow field in different 
regions of the engine. 

During the past two decades, tremendous progress has been made in the field 
of radiative energy transfer in nonhomogeneous nongray gaseous systems. There is 
a renewed interest in investigating various aspects of radiative energy transfer in 
participating medium. Radiative interactions have become important in many engineering 
problems involving high temperature gases. Recent interest lies in the areas of design 
of high pressure combustion chambers and high enthalpy nozzles, entry and rentry 
phenomena, hypersonic propulsion, and defense oriented research. 

Basic formulations on radiative energy transfer in participating media are available 
in standard references [1-5]*. The review articles presented in [6-15] are useful 
in understanding the radiative properties of participating species and the nature of 
nongray radiation. The validity of radiative transfer analyses depends upon the accuracy 
with which absorption-emission and scattering characteristics of participating species 

* The numbers in brackets indicate references 
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are modeled. There are several models available to present the absorption-emission 
characteristics of molecular species and these are reviewed in [12,13]. These models 
have been used to investigate radiative interactions in several duct flows [16-29], 

In a hypersonic propulsion system, the temperature ranges from 1000 - 5000 K. 
In this range various nonsymmetric molecules such as H 2 O, CO 2 and OH become 
highly radiative participating. Infrared absorption and emission of thermal radiation 
is a consequence of coupled vibrational and rotational energy transitions. A diatomic 
molecule is the simplest molecule that will undergo such transition. However, symmetric 
diatomic molecules, such as H 2 , O 2 , and N 2 have no permanent dipole moment and thus 
are transparent to infrared radiation. For unsymmetric diatomic and triatomic molecules, 
such as OH, CO, CO 2 , and H 2 O, the infrared spectrum will consist of fundamental 
vibration- rotation bands occuring at the fundamental frequencies of vibration of the 
molecule, followed by overtone and combination bands [1]. 

In the combustion temperature range, some diatomic and triatomic molecules are 
highly radiative participating species. Various investigators have studied the effect of 
radiative transfer for channel flows. Marring and Hwuang [30] solved the energy equation 
for steam flowing between two parallel black walls. The flow was assumed to be steady 
and the radiation transfer in the flow direction was neglected. It was shown that radiative 
flux peaks at a small distance from the wall, instead of at the wall. This effect was also 
noted by Viskanta in a gray analysis for radiation and convection interaction between 
plates of constant temperatures [31]. This is because at the wall the effect of positive 
radiant heat flux from the wall is partially cancelled by the negative flux from the layers 
of hot gas next to the wall. At a small distance into the stream, however, the flux from 
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the wall and hot gas combine to give a maximum heat flux. 
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Kobiyama et al. [32] studied the problem of combined radiation and convection 
for compressible laminar flow between two isothermal parallel plates. A comparison 
between temperature profiles calculated with the treatment of one and two dimensional 
radiation shows a considerable temperature difference at the enterance region of the 
heating zone. The problem of combined convection and radiation in a rectangular duct 
was studied also by Im amd Ahluwalia [33] for compressible turbulent flows. The 
moment method was employed in this study to solve the radiative flux equation. It was 
concluded that radiative heat transfer causes the thermal boundary layer to grow and skin 
friction to decrease. The velocity profile was not effected by the radiative heat transfer. 

The problem of radiative interaction of gray and nongray absorbing - emitting 
species between two parallel plates and in a circular tube was studied by Tiwari et al. 
[34], In this study, both laminar fully developed incompressible as well as entrance 
region subsonic flows were considered. Results showed that radiative interactions were 
quite significant in fully developed incompressible flows. For subsonic flows, the flow 
field was found to change significantly due to radiative interaction. Tiwari and Singh [20] 
investigated the transient radiative interactions of nongray absorbing emitting species 
in laminar fully developed flows between two parallel plates. The particular species 
considered were OH, CO, CO 2 , H 2 O, and different mixtures of these species. Their 
results demonstrated that H 2 O is a highly radiation participating species compared to 
CO 2 , CO, and OH. The effect of radiation increases with increasing plate spacing, and 
the radiative heat transfer was more pronounced at higher wall temperature and pressure. 
It was also shown that optically thin limits overestimates the influence of radiation. 
Soufiani and Taine [35] studied the H 2 O - air mixture for the above geometry and 
reached the same conclusion . 
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For numerical investigation of chemically reacting flows, an appropriate chemistry 
model must also be selected. Depending on the ratio of the chemical and fluid dynamic 
time scales, the suitable chemistry model could be a frozen flow model, a finite rate 
model, an equilbrium (a complete reaction) model. In general, the finite rate model is the 
most accurate one. In the past several years, a number of finite rate chemistry models for 
hydrogen - air systems have been introduced in the literature. Rogers and Schenayder 
[36] propsed as many as 60 reaction paths in their model; this is certainly one of the 
most complete representatives of hydrogen - air reaction. Unfortunately, the enormous 
number of reaction paths and chemical species involved in the model makes it unfeasible 
for numerical investigation of engineering problems. Intermediate level models are 
reduced to 12 species and 25 reaction paths, and eight species and eight reaction paths. 
Except for some inaccuracies during the ignition delay period, the eight reaction path 
model performs as well as the 25 reaction path model. Although these models are less 
tedious than the 60 path model, they are expected to be too costly for use in routine 
parametric studies. The global two step model proposed by Rogers and Chinitz [37] is 
an inexpensive and attractive model for preliminary investigation of reacting flows. This 
model was deduced by fitting the temperature history of a 28 reaction model [38] used in 
a series constant pressure stream — tube calculations. There are a number of limitations 
to this model, such as ignition phase inaccuracy and a tendency to overpredict the flame 
temperature. This model is considered primarily for initial parametric study. 

A more realistic chemistry model was used by Drummond [23] in numerical 
simulation of a supersonic reacting mixing layer. To explore the behavior of such flows, 
detailed physical models of convective and diffusive mixing and finite rate chemical 
reactions in supersonic flow were developed. The finite rate chemistry model consisted 
of 18 reaction paths and nine species. In this study, two numerical algorithms were 
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constructed to solve the governing equations. The first algorithm was developed by 
modifying the unsplit finite difference scheme of MacCormack. The second algorithm 
employed a hybrid pseudo - spectral technique in the normal to the flow direction for 
improved resolution of the reacting flowfield. The finite difference algorithm was used 
in the streamwise direction. It is important to note that the use of a more complete 
chemistry model rather than the global model in the fluid dynamics equation did not 
result in a set of temporally stiff equations. Incorporation of the finite rate chemistry 
model into the fluid dynamics equation can create a set of stiff differential equations. 
The stiffness is due to a disparity in the time scale of the governing equations. In a 
time - accurate solution, after the large initial transients have decayed and the solution 
is changing slowly, taking a larger time step is necessary for efficiency purposes; but 
explicit methods still require small time steps to maintain stability. An eigenvalue 
problem with stiff ordinary differential equation ( ODE ) has been solved to express this 
clearly in [39]. One way around the problem is to use a fully implicit method. This 
method however requires the inversion of a large coupled system of linear algeabaric 
equations. The use of a semi - implicit technique, suggested by several investigators 
[40] provides an alternative to the above problem. In this technique the source term, 
which is the cause of stiffness, is treated implicitly, and the other terms in the governing 
equations are treated explicitly. 

A comparison of different chemistry models used in investigating radiative 
interactions was done by Chandrasekhar et al. [28]. The results indicated that the 
2 step model predicts ignition before the shock ( shorter ignition delay ) due to high 
temperatures in the boundary layer, while the 18 - step model predicts a longer ignition 
delay. The ignition delay predicted by the 35 - step model appeared to be the average 
of the above two models. A grid sensitivity analysis was also conducted, and it was 
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found that a 31 X 31 grid was the most appropriate grid. 

The purpose of this study is to investigate the effect of radiative heat transfer in 
chemically reacting supersonic flow in various ducts under different physical and flow 
conditions. This provides essential information for investigating the effect of radiative 
interactions in the combustor of a supersonic combustion ramjet (scramjet) engine. 
Several codes have been developed to compute the flow field in a scramjet engine 
[21-24], The combustion of hydrogen and air results in absorbing - emitting gases such 
as H 2 O, OH and NO. Specific attention, therefore, is directed towards investigating 
the radiative contribution of these gases under realistic conditions. Based on the study 
of Ref. 28, a 31 X 31 grid was selected for this study with the 18 step chemistry 
model. The basic code is converted for the axisymmetric case and several parametric 
studies are conducted. 

A brief discussion of the scramjet engine and the physical problem considered is 
presented in Chap. 2 along with the governing equations. Chapter 3 provides the 
formulation of the absorption coefficient and the pseudo-gray model for the absorbing- 
emitting species. The grid generation technique and solution procedure for the governing 
equations are presented in Chap. 4. Discussion of the results for several specific cases 
are provided in Chap. 5, and specific conclusions and recommendations for future 
studies are provided in Chap. 6. 
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GENERAL FORMULATION 
2.1 Physical Model 

As mentioned in the introduction, the scramjet engine has been a candidate for 
propelling hypersonic vehicles. In Fig. 2.1, various air breathing and rocket propulsion 
alternatives are shown. For Mach numbers zero to three, turbojet air-breathing systems 
have the highest performance. Above a free stream Mach number of three, turbine inlet 
temperatures constrain performance, and the ramjet becomes more attractive. At about 
a Mach number of six, the performance of ramjet is greatly reduced. This is due to 
the dissociation of the reaction products, which is caused by slowing the supersonic 
flow to subsonic flow through the normal shock that exist in a ramjet. Therefore it is 
more efficient to allow the engine internal flow to remain at a supersonic speed. Thus 
for Mach numbers of six and beyond, the fixed geometry scramjet is clearly superior 
for propelling a vehicle at hypersonic speeds. Hydrogen has been selected as the fuel 
for the scramjet due to its capability of cooling the engine and the airframe and also 
because of its high impulse level. 

The scramjet engine has several identical modules as shown in Fig. 2.2. The 
forebody of the aircraft acts as an inlet for precompression and the afterbody as a nozzle 
for post expansion. The inlet region starts with the forebody of the vehicle and ends 
with the minimum cross sectional area of each module. In the first part, the air is 
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compressed by the oblique shock generated from the forebody before it enters the engine. 
For numerical solution, the flow is best represented by the Navier-Stokes equations in 
the inlet area of the engine. Using the Euler equations away from the wall region and 
the boundary layer equation near the wall region can be complicated by the oblique 
shock interaction with the boundary layer. This can cause flow seperation and as a 
result flow cannot be represented by these equations. Three-dimensional Navier-Stokes 
equations have been successfully employed by Kumar [21] to investigate the flow field 
in this region. Chitsomboom et al. [24] have employed the parabolized Navier-Stokes 
equations to solve this problem, with limited success. 

The combustor region is the most complex part of the scramjet engine flow field. 
As a result, a great deal of research is directed toward better understanding of the 
combustor flow field. The flow in this region is mostly supersonic, but does have a 
subsonic region near the fuel injection region. The fluid dynamics become complicated 
by the fuel injection, flame holding, chemistry, radiation and turbulence. The flow field 
in this region is represented by the Navier-Stokes equations, turbulence, chemistry and 
radiation modelling equations. Downstream of the fuel injection strut, the flow can be 
represented by the parabolized Navier-Stokes equations [24] 

The physical problems considered here to investigate the effects of radiative 
interactions in supersonic flow are two-dimensional laminar flow between parallel plates 
(Fig. 2.3) and within a circular tube (Fig. 2.4). Another geometry is considered to study 
the effect of shocks and chemical reactions on the radiative heat transfer and this consists 
of a channel with a compression-expansion ramp (Fig. 2.5). The governing equations and 
boundary conditions are provided here for all physical problems considered in this study. 
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(b) Plane radiating layer between parallel boundaries 
Fig. 2.3 Physical model for absorbing-emitting medium. 
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Fig. 2.5 Radiating gas flow in a channel with a compression- 
expansion ramp. 
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2.2 Basic Governing Equation 


The physical problem considered here for basic understanding of radiative interaction 
in supersonic flows is a two-dimensional variable property laminar flow between two 
parallel plates. For this model, the two-dimensional Navier-Stokes equations in fully 
conservative form are used to describe the flow field. These equations, in cartesian 
coordinates, can be written as [25,29] 


dU dF dG 
dt + dx + dy + 

where U, F, G and H are expressed as 


0 


( 2 - 1 ) 
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The viscous stress tensor in the F and G terms of Eq. (2-1) are given as 


( du d\ 

Txy = h W + 


(2-2b) 


f d u dx\ dv 

Tyy - -P + A^ + -^)+2n— (2-2C) 

The quantities q cx and q cy in F and G terms are the components of conduction heat 
flux and are expressed as 


Qcx — 
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T 
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It should be noted that D represents the binary diffusion coefficient and is used for 
all species. Assuming that the Lewis number ( a / D ) is unity, Eq. (2-3) reduces 
to (see Appendix A) 
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(2-4a) 
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7 fi de 

pr dy 


(2-4b) 


where e = h - P Ip. The molecular viscosity, /z e , is assumed to be only temperature 
dependent, and is evaluated from the Sutherlands formula as 


( T T 0 + S 

/ie - M \^ T( J T + g 


(2-5) 


where /j.q and Tq are reference values for individual species and S is the Sutherland 
constant. Here the reference values were selected for pure air because the flow is 
dominated by nitrogen. The total internal energy E is given by 



u z + v * 
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m 


+ £ hifi 

i=l 


( 2 - 6 ) 


Equation (2-1) can be used to obtain solutions for various kinds of compressible 
flows. However, boundary conditions and numerical procedures for different flows are 
quite different. For supersonic flows, the inflow conditions are specified and outflow 
conditions are obtained by extrapolation. The boundary conditions used along solid 
surfaces are u=0, v=0, ^ = 0, and T = T w . 

The governing equations and boundary conditions for the supersonic flow through a 
channel with a compression-expansion ramp are essentially the same as for the parallel 
plate geometry. However, a strong shock is produced at the compression comer and the 
flow becomes highly reacting from the begining of the X 2 -coordinate (Fig. 2.5). 
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The basic governing equation for chemically reacting compressible flow through a 


circular tube can be written as [39] 


OT 1 a (yg) H 
dt dx y dy y 


(2-7) 


where x, y represent streamwise and radial coordinates, respectively. In Fig. 2.4, the 
radial coordinate is denoted by r. The definitions of U, F and G in Eq. (2-7), are same 
as given in Eq. (2-1) but vector H is expressed as 
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0 
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P + A 06 

0 


( 2 - 8 ) 


The viscous terms in this case are given by 
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(2-9d) 


. f du du v\ v 

m = H^ + a? + yJ + "y 

The boundary conditions for the circular geometry are similar to the parallel plate 
geometry. 


2.3 Chemistry Model 


In the present work, the finite rate chemical reaction of gaseous hydrogen fuel and 
air is studied. The reaction is modeled by a nine species and eighteen reaction model 
presented in Table 2.1. Eight of the chemical species ( H 2 , O 2 , H 2 O, OH, H, O, HO 2 , 
h 2 o 2 ) are active and the ninth ( N 2 ) is assumed inert. The forward reaction rate of 
each reaction j is given by the Arrhenius law 



( 2 - 10 ) 


Values for A, N and E are given in Table 2.1. Knowing the forward rate constants and 
using the equilibrium constant, the backward rate constants can be obtained from 


k 


bj 
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The equilibrium constant defined in the abve equation is given by 
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Once the forward and reverse reaction rates have been determined, the production 
rates of eight species can be found from the law of mass action. For the general 
chemical reaction 


N, N. 

^7ji c i <=> ^7j"i c i , j = 1,2,3, N r (2-13) 

i=l i=l 

the law of mass action states that the rate of change of concentration of species i 
by reaction j is given by 



N s , N, „ 

k i. II <V - II C> 

i=l i=l 


i = l,2,...N s (2-14) 


The net rate of change in concentration of species i by reaction j is then found by 
summing the contribution from each reaction as 


Ci = Y, < C i)j (2-15) 

i=l 

Finally, the production term, w, present in the governing equation is obtained from 

Wi = CjMi ( 2 - 16 ) 

where Mj is the molecular weight of the i^ 1 species. 

For some reacting flows, the Jacobian of the chemical system must be evaluated 
with respect to the field variables, (including each species). Certain chemical systems 
are numerically unstable when integrated with a time step comparable with the CFL 
condition. One remedy for this instability is treating implicitly the chemical source terms 
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Wj. This technique requires the Newton linearization of the chemical source term about 
the previous time step and thus results in the need for evaluation of the Jacobian matrix 
which relates the rate of change of each source term with respect to each species. This 
Jacobian can be calculated either numerically or analytically, and presents no special 
problem once a system of reactions and species has been introduced. 

Analytically, the chemical Jacobian is obtained by differentiating Eq. (2-14), and 
is defined by 


5Ci 

<9C m 


N s 


E (*$ - ?ii) 


KfV N > , K 


b ‘ 7 i™ jp c j> 


J=1 


i=l 


i=l 


i = 1,N S 

(2-17) 


The numerical Jacobian is calculated by approximating the limit formula for the 
derivative as 


dCj = lim Cj (P, T, Cl, C m + AC m , C m+1 , ..C ns ) 

dC m AC..,— o AC m 
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Table 2.1. Hydrogen-Air Combustion Mechanism [22] 


— 


REACTION 

A(moles) 

N(cm 3 ) 

E(calories/gm-moIe) 



** following reactions constitute the 18-step model ** 



0) 

0 2 + H 2 — OH + OH 

1.70 x 10 13 

0 

48150 


(2) 

0 2 + H — OH + 0 

1.42 x 10 14 

0 

16400 

— 

(3) 

H 2 + OH — H 2 0 + H 

3.16 x 10 7 

1.8 

3030 


(4) 

H 2 + O — OH + H 

2.07 x 10 14 

0 

13750 


(5) 

OH + OH — H 2 0 + O 

5.50 x 10 13 

0 

7000 


(6) 

H + OH + M — H 2 0 + M 

2.21 x 10 22 

—2.0 

0 


(7) 

h + h + m->h 2 + m 

6.53 x 10 17 

—1.0 

0 


(8) 

H + 0 2 + M — H0 2 + M 

3.20 x 10” 

—1.0 

0 

: : 

(9) 

OH + H0 2 - 0 2 + H 2 0 

5.00 x 10 13 

0 

1000 


(10) 

H + H0 2 — H 2 + 0 2 

2.53 x 10 13 

0 

700 


(11) 

H + H0 2 — OH + OH 

1.99 x 10 14 

0 

1800 


(12) 

O + H0 2 - O z + OH 

5.00 x 10 13 

0 

1000 

«pii 

(13) 

H0 2 + H0 2 - 0 2 + H 2 0 2 

1.99 x 10 12 

0 

0 


(14) 

h 2 + HO z — H + H 2 0 2 

3.01 x 10 u 

0 

18700 

i : ? 

(15) 

OH + H 2 0 2 — H 2 0 + H0 2 

1.02 x 10 13 

0 

1900 


(16) 

H + H 2 0 2 — H 2 0 + OH 

5.00 x 10 14 

0 

10000 


(17) 

o + H 2 0 2 — OH + H0 2 

1.99 x 10 13 

0 

5900 


(18) 

h 2 o 2 + m — oh + oh + m 

1.21 x 10 17 

0 

45500 


** remaining reactions complete the 35-step model ** 


— 

(19) 0 2 + M — O + O+ M 

2.75 x 10 19 

—1.0 

118700 


(20) N 2 + M-*N + N + M 

3.70 x 10 21 

—1.6 

225000 


(21) N + O 2 — * 0 + NO 

6.40 x 10 9 

1.0 

6300 

* j 

(22) N + NO — O + N 2 

1.60 x 10 13 

0 

0 


(23) N + OH — H + NO 

6.30 x 10" 

0.5 

0 


(24) H + NO + M — * HNO + M 

5.40 x 10 15 

0 

—600 


(25) H + HNO — H 2 + NO 

4.80 x 10 12 

0 

0 


(26) O + HNO — OH + NO 

5.00 x 10 u 

0.5 

0 

— 

(27) OH + HNO — H 2 0 + NO 

3.60 x 10 13 

0 

0 


(28) HOi + HNO — H 2 Oi + NO 

2.00 x 10 12 

0 

0 

. , 

(29) H0 2 + NO — OH + N0 2 

3.43 x 10 12 

0 

—260 

K :S 

(30) H + NOi — OH + NO 

3.50 x 10 14 

0 

1500 


(31) O + N0 2 — 0 2 + NO 

1.00 x 10 13 

0 

600 


(32) N0 2 + M — O + NO + M 

1.16 x 10 16 

0 

66000 


(33) M + OH + NO — HN0 2 + M 

5.60 x 10 1S 

0 

—1700 


(34) M + OH + NOj — HN0 3 + M 

3.00 x 10” 

0 

—3800 


(35) OH + HN0 2 - H 2 0 + NCh 

1.60 x 10 12 

0 

0 

— 

** following reactions constitute the global 

2-step model [4, 16, 

18, 19] ** 


(1") H 2 + 0 2 — 2 OH 

11.4 x 10 47 

—10.0 

4865 

— 

(2") 2 OH + H 2 — 2 H 2 0 

2.50 x 10 64 

—13.0 

42500 
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2.4 Thermodynamic model 


The specific heat of individual species, Cpj , is assumed to be a linear function 
of temperature, i.e., 

C Pj = ajT + bj (2-19) 

where aj and bj are constants which are obtained by curve fitting the thermodynamic 
data of Ref. 40. The numerical values of these constants are given in Table 2.2. The 
specific heat of the mixture is computed by summing specific heats of individual species 
weighted by species mass fraction 


m 

cp = E °pj f j 

j=i 

The static enthalpy of the mixture can be expressed as 


m 


>? + / Cpj dT 
To 

The total enthalpy can now be evaluated as 


E 

j=l 


( 2 - 20 ) 


( 2 - 21 ) 


H = h + 0.5(u 2 +v 2 ) 

Combining Eqs. (2-19), (2-21) and (2-22) the total enthalpy is expressed as 


( 2 - 22 ) 


H = E 


a*T 2 

h? + + b;T 


j=l 


fj + 0.5 (u 2 + v 2 ) 


(2-23) 


where h°j is the sensible enthalpy of individual species at a reference temperature Tq 
= 0 K . The gas constant for the mixture is evaluated by a mass weighted summation 
over all of the species as 
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( 2 - 24 ) 


R = £ f j R j 

j=i 

The equation of state for the mixture of the gases therefore can be written as 


P = pRT 


( 2 - 25 ) 
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Table 2. 2 Numerical values of various constants 


Species H(J/Kg) a b 


0 2 

-271267.025 

0.119845 

947.937 

h 2 o 

-13972530.24 

0.43116 

1857.904 

h 2 

-4200188.095 

2.05960 

12867.46 

OH 

1772591.157 

0.16564 

1672.813 

N 2 

-309483.980 

0.10354 

1048.389 

NO 

2733268.38 

0.34206 

1256.78 
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Chapter 3 


RADIATION TRANSPORT MODEL 

3.1 Radiation Absorption Model 

In order to include the effects of radiative interaction in a physical problem, it is 
essential to accurately model the absorption - emission characteristics of participating 
species and provide a correct radiative transfer formulation for the physical process. 

Many models are available in the literature ; to present the absorption - emission 
characteristic of molecular species; a review of important models are available in [13]. 
Perhaps the simplest model is the gray gas model, where the absorption coefficient 
is assumed to be independent of the wave length. Many nongray models are also 
available in literature. In this study the gray gas model is utilized primarily because 
of its simplicity. 

The total band absorptance of a vibration - rotation band is given by 

OO 

A = j [1 - exp(/c w X)d(u; -wo)] (3-1) 

— OO 

where k w is the volumetric absorption coefficient, w is the wave number, wo is the 
wave number at the band center, X = Py is the pressure path length, and the limits of 
integration are over the entire band pass. Various models are used to obtain the relation 
for A in Eq. (3-1). A convienent model to represent the average coefficient of a gray 


25 


gas is the Planck mean absorption coefficient, /c p , which is defined as [1] 

OO 

K P = J K u e bu> (T) du>/ e b (T) (3-2) 

0 

For a multiband gaseous system, this is expressed as 


K p 



£«ui(T)S;(T) 


(3-3) 


where Pj is the partial pressure of j* species in a gas mixture, e wi ( T ) is the Planck 
function evaluated at the band center, and Sj( T ) is the integrated band intensity 
of the i^band. 

As defined in Eq. (3-2), k p is a property of the medium. When « p is evaluated 
at the temperature of the gas, it is usually a mean emission coefficient and it becomes 
equal to the actual mean absorption coefficient only for the conditions of an equilibrium 
radiation field. For a nonuniform temperature field, the mean absorption coefficient 
which is used for the optically thin radiation is the modified Planck mean absorption 
coefficient, which for black bounding surfaces is defined as 


K m (T ; T w ) — 


J (T) (T w ) du> 


(3-4) 


e b (^w) 

Note that «; m is a function of both the gas temperature and the wall temperature. An 
approximate relation between k p and K m is available for infrared radiation as [1,9] 


K m (T, T m ) 



(3-5) 


This expression is usually employed in gray gas radiative energy transfer analyses. 
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Several models for the mean absorption coefficient are available in the literature 
[1,31]. Since these models account for detailed spectral information of molecular bands, 
this approach of radiative formulation is referred to as the “ pseudo - gray formulation”. 
The gray gas formulation for radiative transport is very useful in parametric studies. 

There are several nongray models available in the literature to represent the absorption 
- emission characteristics of vibration - rotation bands. These are classified typically in 
four types: ( 1 ) line - by - line ( LBL ) models, ( 2 ) narrow band models, ( 3 ) wide 
band models, and ( 4 ) band model correlations. A complete discussion on usefulness and 
application of these models is provided in [12,13]. For many engineering applications, 
wide band model correlations provide quite accurate results. The most widely used wide 
band models are those of Edwards [5,1 1] and Tien and Lowder [9]. The Tien and Lowder 
correlation for the total band absorptance is a continuous correlation and is given by 


A (u, /?) = A ( u , /?) / A 0 = ln{uf(t)[(u + 2)/(u + 2f(t))] + 1} (3-6) 

where 

f (t) = 2.94 [1 - exp (-2.60t)] , t = 0/2 (3-7) 

and u = SX/Aq is the nondimensional path length, /? = 27T7/d is the line structure 
parameter, 7 is the line half width, S is the integrated band intensity, and Aq is the 
band width parameter. Equation (3-6) provides accurate results for pressures higher 
than 0.5 atmosphere [12,13]. 

Spectral properties and correlation quantities for various radiation participating 
species are available in [5,9,11]. These are useful in gray as well as nongray radiative 
formulations. 
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3.2 Evaluation of Planck Function and 
Planck Mean Absorption Coefficient 

The spectral emissive power of a black body at a temperature T is given by the 
Planck’s law as [1] 

27rhi/ 3 n 2 


e b (^ T) 


L 0 


ex P (gp) - 1 


(3-8) 


where 


co = speed of light in vaccum (2. 998 x 10 10 cm/sj 
h = Planck's constant ^6.625 x 10 _27 erg — sj 
k = Boltzman constant (l.380 x 10 -16 erg/Kj 
n = Index of refraction of the medium 


(3-9) 


v = Frequency of radiation 

The speed of light c in a medium is expressed as c = A u, where A is the wavelength; 
this is related to the speed of light in a vaccum by c = cg/n. Thus, v = co/(nA). If it 
is assumed that n is independent of frequency, then 


Since e^dA = — e^di/ , Eq. (3.8) can be expressed as 

Ci 


where 


e b (A,T) = 


n 2 A 5 


exp 


(to- 


Ci = 27ihC 2 = ^374.135577 x 10 _7 erg - cm 2 /s) 


(3-10) 


C 2 = hC 0 /k = (l43. 9257246 x 10 _2 cm - k) 

The index of refraction n is assumed to be independent of frequency and its value 

is taken to be unity. Consequently, Eq. (3-10) is applicable to most gaseous systems. 
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It is often convenient to express Eq. (3-10) in terms of the wave number u;=l/A. 
In order to accomplish this, one should note that dw=-A _2 dA and e b ^ dA=-e bu ,dw, such 
that %w = /\ 2e bA- Thus Eq. (3-10) can be expressed as 


e b( w > T ) 


27rhc 2 cj 3 
exp (Iff) - 1 


Cjo ; 3 

exp (C2^) — 1 


(3-11) 


From the definitions of Ci and C 2 , and u expressed in cm l , it is evident that the units 
for the Planck function in Eq. (3-11) are erg/s-cm or W/cm. 


For certain applications, it is desirable to have a convenient expression for the 
temperature derivative of the Planck function. The quantity of frequent interest is 
de b (t*,’,T)/dT and this can be obtained from Eq. (3-11) as 


de b (u;,T) 

dT 


cic 2 ^ 4 exp (^) 

T 2 [exp (^) - l] 2 


( 3 - 12 ) 


The Planck mean absorption coefficient k p is defined by Eq. (3-2) and is expressed 
here as 


ftp 


OC 

/ K w e b ( w iT) du> 
0 


aT 4 


(3-13) 


It should be noted that for an isotropic medium k u depends only on temperature; i.e., 
k w =«: u ,(T). Thus, Kp is a property of the medium, provided that conditions of local 
thermodynamic equilibrium exist [1]. 

The integrated intensity of a band S (i.e., the integrated band absorption) is defined as 
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(3-14) 


S(T)= j (^)du. 

A a; 

where Aw represents the spectral range of the band. If the value of S is available 
at a given (standard) temperature Tq, then the value of S at any other temperature is 
obtained from the relation 


S(T) = 0r)s(T o )F(T) (3-15) 

where F(T) is a function of temperature and its value depends upon the nature of the 
band. For fundamental vibration-rotation bands, the value of F(T) is unity, i.e., F(T)=1. 
For combination and overtone bands of certain species, the values of F(T) are available 
in the literature. However, for many species, the relationship for variation of S with 
temperature for combination and overtone bands are not yet available. 

Within a narrow vibration-rotation band, the Planck function variation with the 
wavelength is small and its value may be evaluated at the band center, i.e., eb( u>, T ) = 
e^( w c , T ). With this assumption, Eq. (3-13) is written for a multiband system as 


«p(T) = 


N 

E 

i=l 


e b(^ T ) f 


(3-16) 


p ^' a T 4 

where N represents the number of vibration - rotation bands of a gas. By employing the 
definition of S given in Eq. (3-14), Eq. (3-16) is expressed for a single homogeneous 
gas as 


P E [e b (w c ;, T) S; (T)] 

«p(T) = —4 (3-17) 
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where Si( T ) represents the integrated band intensity of the i th band of the gas. Equation 
(3-17) can be modified to apply to a mixture of different gases as 

E p j j.E [ e b Ki. T) Si (T)] | 

«P (T) = - gT4 1 (3-18) 

where j denotes the number of species in the mixture and Pj is the partial pressure 
of the j^species. 


3.3 Radiative Flux Equation 

The equations of radiative transport are expressed generally in integro - differential 
form; the integration involves both the frequency spectrum and the physical coordinates. 
To overcome the complexity of the radiative transport equations, a tangent slab 
approximation was employed in [26,27], This approximation treats the gas layer as a 
one dimensional slab in the evaluation of the radiative flux. The multi - dimensional 
equations of radiative transfer are formulated for an arbitrary geometry, and then an 
approximate method is used to present the formulation for gray and nongray gases. 
The complete formulation is available in [27], and essential relationship are provided 
in this section. 

3.3.1 Basic Formulation 

The radiative transport equations in the present study are obtained only for an 
absorbing - emitting medium contained within solid walls of different configuration as 
shown in Figs. 2. 3-2.5. The general formulation of radiative transfer for the gas under 
the condition of local thermodynamic equilibrium is given by (see Figs. 2.3 and 3.1) 
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Fig. 3.1 Coordinate system for 


one-dimensional radiative transfer. 
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(3-19) 


/ MOd* 7 - 7 MOd* 

Ii/(P) = Ii/(Pw) e° + J (P') I b „ (P') e « dr 7 

0 

The first term on the right hand side represents the contribution from the wall to the 
intensity at P, and the second term represents the contribution from the intervening gases 
between P and P w - The origin of the coordinate system in Eq. (3-19) is chosen at the 
point P. The radiative flux at the point P in the direction of l is 


oo 


qr 


-IS 


I dI dv dfi 


(3-20) 


47r 0 


Substituting Eq. (3-19) into Eq. (3-20), the radiative flux is expressed as 


qR 


t °f -/«.(*>< ^ 

/ /l„(P„)e • IdvASl 


47 r 0 
oo r w 


(3-21) 


-f. 


+ IJ I K V (P 7 ) I b „ (P # ) e " fdr 7 dv dfi 


4ir 0 0 

The divergence of the radiative flux is formulated as 


v »qR 


oo 


II 

47 r 0 


l . VI„ dv dn 


(3-22) 


A beam of intensity I^( r, Q ) travelling in the l direction satisfies the equation of 
radiative heat transfer 


/ . VI*, = 7 „I„ + K v I hv ~ Wv (3-23) 

where is the scattering coefficients and /3 U = 71 , + is referred to as the extinction 
coefficient. For negligible scattering ( ju = 0 ), Eq. (3-23) reduces to 
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(3-24) 


l ■ VIj / — Ki/ (I|jj/ ~ Ij /) 

By combining Eq. (3-22) and Eq. (3-24), one obtains 


Now substituting Eq. 
expressed as 


V • q R = 

(3-19) into 


OO 

1 1 — ^*0 dz/ (3—25) 

47 r 0 

Eq. (3-25), the divergence of the radiative flux is 


oo oo 

v qR. = 47T J k u (p) I bl/ (p) du — J J K l/ (p)I l/ (p Vf ) 


f 


e 0 


4t r 0 


(3-26) 


OO 


f f f 

/ / / (p) «j/ (p ) Ibt/ (p') e 0 dr 7 du dfi 


4 tt 0 0 

Equations (3-21) and (3-26) are used to obtain various approximate forms for the 
radiative flux and its divergence. 


3.3.2 Gray Formulation 

In the previous section, it was observed that the radiative flux terms are represented 
by an integral equations. Solving these equations is extremely time consuming because 
of the complexity of integration over space and frequency. Therefore, a pesudo - gray 
model is selected for efficient parametric studies. To express the radiative flux for a 
gray medium, one may assume that is independent of the frequency. This is rarely 
a physically realistic approximation; but it serves as an initial stepping stone towards 
nongray analyses. Therefore, Eq. (3-21) and Eq. (3-26) for a gray medium are written as 
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/ «(0d* 


/ - j r r 

I(P w )e • df*+ / / k (p 7 ) I b (P 7 ) e « dr'dn (3-27) 
4?r 4 tt 0 


- f «({)* 

e u 


and 


r ~ f «(O d £ 

V • q R = 4*7e(p)I b (p) - j «(p)I(P w )e 0 dfi 


4tt 


r w 

~ 1 1 K(p ) K ( p/ ) Ib ( p, ) e 
4?r 0 


/ *(O d £ 

« dr' dft 


(3-28) 


To solve the radiative flux terms for the gray medium, Eqs. (3-27) and (3-28) can 
be transformed into Cartesian coordinates, and then any one of the standard integration 
techniques to evaluate the radiative flux term are applied. It is more convenient and 
efficient to convert the equations to a set of ordinary differential equations (ODE). For 
the present case, differentiating Eq. (3-28) by using the Leibnitz rule results in 


V 2 q R = 4K(p)^ + 


J * 2 (p)I(P„) 

4tt 


e 


f «(0d* 

® dO 


+ J / « (P) « (p') r b 

4t r 0 


r' 

- / 

(p 7 ) e o dr 7 d Q 


(3-29) 


A substitution of Eq. (3-29) into Eq. (3-27) gives a second order nonhomogenous 
ordinary differential equation for one dimensional radiative transfer as 
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(3-30) 


d2 qR 

dr 2 


« 2 (P)qR 


4«(P) 


de b 

dr 


It should be pointed out that if the tangent slab approximation is employed, then the 
method of exponential kernal approximation is used to convert the equations to a set 
of ODE’s. The transformation of Eq. (3-28) into a differential equation by exponential 
kernal approximation is shown in Appendix C. The result is of the same form as Eq. 
(3-30), but the coefficients are different; i.e., 


d 2 qR 

dr 2 


j K 2 (P)q R + Sk(P)^ 


(3-31) 


Equations (3-30) and (3-31) differ in the coefficients due to the employment of the 
exponential kernal approximation in the second approach. These equations require 
two boundary conditions. For nonblack diffuse surfaces, the boundary conditions 
corresponding to Eqs. (3-30) and (3-31) are given in [12]. 

For black walls and for Tj = T 2 , the boundary conditions for Eq. (3-31) become 



2qR(°) 


J_ f d qR ^| 

T 0 V d£ /£ =0 


(3-32) 


where 


to — KpL , 



(3-35) 


For a black circular tube, the spectral radiative heat flux in the radial direction is 
given by the expression [17] 
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qR(r) 


jr 

2 r 


7 j { J F w (r') k^aexp 


r SU17 


bkg, 
cos 7 


(r - r') 


di^ 


Tu 

J F w (r') k^aexp 

r 


bkg; 

COS 7 


(r - r) 


dr' 


(3-33) 


U 

+ / F w (r') k^aexp — (r + r' — 2rsin7) dr' 

J L cos 7 


}d7 


r sin7 


where Fu, (r ) = e w (r ) - eu,(T w ) and constants a and b have values of unity and 5/4, 
respectively. A combination of Eq. (3-20) and (3-34) provides a convenient expression 
for the total radiative flux for nongray analyses. 

For a gray medium, the expression for the total radiative flux can be obtained for 
a circular tube from differential approximation as [1,17] 


d_ 

dr 




- 7 k pqR = 3<Tk i 


dT 4 

dr 


(3-34) 


The boundary condition for this equation is found to be 


?qR(i) = -r- 


737 (£qn) 


to U d £ 


-U=i 


> qR(°) = ° > 


where 


TO — K p r 0 , £ — 


ro 


(3-35) 


(3-39) 


With certain modifications, the radiative flux equations presented in this section can 
also be used to investigate the radiative interactions in the flow direction. The procedure 
for doing this is provided in [27]. 
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Chapter 4 


METHOD OF SOLUTION 

The grid generation technique and solution procedures for the governing equations 
using the unsplit MacCormack [43] technique are briefly described in this chapter 

4.1 Grid Generation 

Grids are generated using an algeabaric grid generation technique developed by 
Smith and Weigel [42], From the computational point of view, it is desirable to have 
a uniform rectangular grid enclosed in a cube, where the exterior of the grid represents 
the physical boundaries. To have such grids, the body - fitted coordinate system is 
transformed linearly from the physical domain ( x, y ) to the computational domain ( £, 
i] ) as follows : 


Lower Boundary 


Xi = X(£,0) 
Y i = Y(£,0) 


(4-1) 


Upper Boundary 


X 2 

Y 2 


X(G1) 

Y (£, 1) 


(4-2) 


Between the Boundaries 


X = X (£, 1) 77 + X (£, 0) (1 — rf) 
Y = Y(Gl)i? + Y (GO) (!-»/) 


(4-3) 
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where 


0 < £ < 1 ; 0 < 77 < 1 

The grid should be concentrated in the region of high gradients to accurately predict 
the solution. Therefore, more grid points are required near the solid boundaries. The 
concentration of the grid in the 77 direction can be accomplished by 


_ (fty + 1) ~ (fty ~ 1) exp [-c (77 - 1 + a) / (1 - a)] 
(2a + 1){1 + exp [c (77 - 1 + a) / (1 - a)]} 


(4-4) 


where 



If a is equal to zero , the compression takes place only near the lower wall ( 77 = 0 


), and if a is set equal to one half, the compression takes place near both the walls. 
The term /3y has a value between one and two, and as it gets closer to one, the grid 
becomes more concentrated near the walls. Employing this concentration, Eq. (4-3) 
is written in terms of 77 as 

X = X(i,l)fj + X(£,0)(1-t7) 


Y = Y (£, 1) 77 + Y (£, 0 ) (1 — fj) 


(4—5) 


where 


0 < 77 < 1 

It should be noted that the grid is concentrated near solid walls in the normal direction 
to resolve the boundary layer. Uniform spacing is used in the streamwise direction. 


4.2 Solution of the Governing Equations 

The MacCormack’s finite-difference scheme is used to discretize and solve the 
governing equations. As such, Eq. (2-1) is expressed in the computational domain as 
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+ H 


0 


( 4 - 6 ) 


where 


dt [ 
dt 


dF dG 

d( drj 


U = UJ 


F = FY,, - GX^ 
G = GX^ - FY^ 


H = HJ 


J = X^Y,j — 

Equation (4-6) is discretized temporally and written as 


U n+1 


U n - At 


dF n 

d( 


The source term H n+1 must be linearized next, 
time to give 


+ 


+ H n+ 1 

or) 


( 4 - 7 ) 


It is expanded in a Taylor series in 


H n+1 = H n + At-^- + 0 (At) 2 


or 


H n+1 


H n + At 


m 

au 


/ u“+i 
At 



(4-8) 


A substitution of Eq. (4-8) into (4-7) gives the temporally discrete equation in delta 
form as 


’ . dH 

I T At — — 

AU n+1 = 

-At 

'dF dG - 
4 — 5 — T H 

< 9 U_ 



d£ dr) 


( 4 - 9 ) 


where U n + 1 — U n is expressed as A U 


n + 1 


3H 


, is the Jacobian matrix of H and 


’ IduJ’ 

[I] is the identity matrix. The terms of the Jacobian matrix are detailed in Appendix D. 
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Once the temporal discretization used to construct Eq. ( 4-9) has been performed, 
the resulting system is spatially differenced using the unsplit MacCormack predictor 
- corrector scheme. This results in a spatially and temporally discrete, simultaneous 
system of equations at each grid point. Each simultaneous system is solved using the 
Householder technique [45] in combination with MacCormack technique, which is then 
used to advance the equation in time. The modified MacCormack scheme then becomes 


1+ At 



Aun +1 = -At 


dF dG 
+ 3^ + H 


1 n 


U 


(4- 10a) 


U” +1 = UR + AUR +1 (4- 10b) 
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(4- 10c) 


UR +1 = UR + 0.5 [auR + 1 + UR +1 ] (4— lOd) 

Equations (4-10) are used to advance the solution from time n to n + 1, and iteration 
process is continued until a desired integration time has been reached . 

The details of gray gas as well as nongray radiative flux formulations and solution 
procedure are available in [1], For the gray gas model, the governing equations are 
discretized using a central difference scheme. The discretization of Eqs. (3—31) and 
(3-34) results respectively in 
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q j+1 = RHS (4-11) 
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e j+i ~ e j , e j ~ e j-i 

^j A yj A yj 


, Ayj = yj - yj-i , Pj 


yj+i ~ yj 
yj-yj-i 


Equations (4-11) and (4-12) along with boundary conditions given by Eqs. (3-32a) 
and (3-32b) form tridiagonal systems of equations which can be solved efficiently by 
the Thomas algorithm. 

In the nongray gas formulation, the divergence of the radiative flux is evaluated 
using a central difference scheme and is treated as radiative source term in the energy 
equation. Since the radiative flux terms involves integral formulation, unlike other flux 
terms which are only in a differential form, it is uncoupled and treated seperately. 

4.3 Physical Conditions and Data Source 

The physical conditions for which specific flowfield analysis and computations are 
needed are discussed in [21-28]. In this work selected parametric studies have been 
conducted for certain flow and physical conditions. Radiation participating species which 
are considered are H 2 O, OH and NO. Radiative properties of these species are available 
in [5,9,11-13]. Different amounts of these gases, in combination with air, are considered 
for parametric studies. Essential data for the chemistry model employed are obtained 
from Refs. 28 and 29 and these are provided in Table 2.1. 
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For basic studies, the physical dimensions considered for the channel are L=3 cm 
and L x =10 cm, and for the circular tube they are L=D=3 cm and L x =10 cm. Reacting 
flow calculations are obtained for a 10° compression nozzle. 

As discussed in Chap. 2, the governing equations Eq. (2-1) requires boundary 
conditions along all four boundaries. For the cases considered in this study, the inflow 
boundary is always supersonic, so the velocity and species are specified there. The upper 
and lower surfaces are solid boundaries, so the no slip boundary conditions (u=0, v=0) 
are used to specify the velocity components. In addition, along these solid boundaries the 
temperature gradient is held constant for one case and for another case the temperature is 
specified. The outflow boundary is also supersonic, and therefore velocitiey components, 
static temperature and pressure, and species concentration are extrapolated. 

The governing equations also require initial conditions. The equations are initialized 
by setting values of the velocity, static temperature, and pressure, and species 
concentration through the domain to the values at the inflow boundary. Having 
specified all required initial and boundary data, the equations are marched in time from 
the initial time to some final specified time level. 
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Chapter 5 


Fi 


RESULTS AND DISCUSSION 

Based on the theory and computational procedure described in the previous sections, 
an existing code was modified to solve the two-dimensional Navier-Stokes equations 
for radiating supersonic laminar flow between two parallel black plates. A similar code 
was developed for radiating supersonic flows in a circular tube. In most cases, the 
radiative interaction was considered only in the normal direction. Extensive results have 
been obtained for pure H 2 O, OH, and NO as homogenous participating species, and for 
different mixture of these species with air. 

As mentioned earlier, the Planck mean absorption coefficient « p (or K m ) is considered 
to be an optically thin radiation absorption coefficient, although it has been used in other 
optical ranges as well [1,9]. The appropriate absorption coefficient for the optically thick 
radaiation is the Rosseland mean absorption coefficient /cr. It has been noted in [1] that 
if the medium is gray, then Kp = kr = k; otherwise Kp > kr. Thus use of /c p (or K m ) in 
pseudo-gray gas formation will provide maximum influence of radiative interaction. The 
Kp values for H 2 O, OH, and NO have been calculated from Eq. (3-17) by employing 
radiative properties available in [5,9,11], and these are illustrated in Figs. 5.1 and 5.2. 
Figure 5.1 shows the results of 100% homogenous species, whereas results in Fig. 5.2 are 
for different mixtures. The results provide indication of the radiative ability of different 
species at a given temperature. Values of k p for other species are available in [1,9]. 
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Fig. 5.1 Variation of Planck mean absorption coefficient for H 2 0, NO and OH 
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Temperature 

Fig. 5.2 Variation of the Planck mean absorption coefficient for different mixtures of 
H 2 O, NO and OH. 


For basic understanding of supersonic entrance region flow, illustrative results for 
flowfield variables and heat transfer have been obtained for nonreacting flows between 
two parallel plates and within a circular tube. The physical dimensions considered for the 
channel are L=3 cm and L x =10 cm. In most cases, the inflow conditions considered are 
P 0 o=l atm, Too=l,700 K, £/oo=2574 m/s (Mqq-3.0) with varying amounts of radiation 
participating species in combination with air. The surface temperature varies with x, and 
is determined at the first grid point normal to the surface. Certain variations in physical 
and inflow conditions are also considered for parametric studies. The chemical reactions 
are not considered in obtaining the results, and radiative flux results are presented only 
for the normal direction. 

Results for the radiative flux as a function of nondimensional position, are illustrated 
for P = 1 atm. The results presented in Fig. 5.3 for different water vapor concentrations 
indicate that radiative interaction increases slowly with an increase in the amount of the 
gas. It is noted that the radiation flux is approximately zero in the center of the channel 
( y = 1.5 cm ) and is significantly higher towards the top and bottom of the plates. This 
however would be expected because of the symmetry of the problem and the relatively 
high temperatures near the boundaries. The results for 50% H 2 O are presented in Fig. 
5.4 for two different pressures (P = 1 and 3 atm) and x-locations (x = 5 and 10 cm). It 
is noted that the increase in pressure has dramatic effects on the radiative interaction. 

For a mixture of 50% H 2 O in air, the conduction and radiation heat transfer results 
are compared in Fig. 5.5, for P = 3 atm and for two different x-locations (x = 5 and 
10cm). The results demonstrate that the conduction heat transfer is restricted to the 
region near the boundaries, and does not change significantly from one x-location to 
another. The radiative interaction, however, is seen to be important everywhere in the 
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Fig. 5.3 Radiative flux vs y at the channel exit = 3.0 
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Fig. 5.5 Radiative and conductive fluxes vs. y for P = 3 atm, 
and 10 cm, 50% H 2 0, = 3.0. 
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channel, and this can have a significant influence on the entire flow field. 

Comparitive results for 100% homogeneous species of H 2 O, NO, and OH are 
illustrated in Fig. 5.6, at the exit plane (x = 10 cm). As would be expected, the radiative 
contribution of H 2 O (with five bands) is significantly higher than NO and OH. Only 
the fundamental bands of NO and OH are considered in this study, and it is noted that 
NO is a better radiating gas in comparison to OH. 

For a mixture of 25% H 2 O in air, radiative flux for two different plate spacings (L 
= 3 and 6 cm) are illustrated in Fig. 5.7 for two x-locations (x = 5 and 10 cm). The 
rate of radiative transfer is a strong function of the amount of the participating species 
and the pressure path length PL. Consequently, the results for the larger plate spacing 
indicate significantly higher radiative interactions, .. 

The effect of increased Mach number on the radiative transfer is illustrated in Fig. 
5.8 for pure H 2 O and for a mixture of 50% H 2 O in air. The results shown are for 
the exit plane (x = 10cm). At higher Mach number, the boundary layer is relatively 
thinner and the temperature in the boundary layer is significantly higher. This in turn 
results in a higher rate of radiative transfer. 

For a mixture of 50% H 2 0 in air, comparitive results for the parallel plate channel 
and the circular tube are presented in Fig. 5.9, for two x-locations (x = 5 and 10 cm). 
The results for the circular tube exhibit the same trend as for the parallel plate geometry. 
Since the circular geometry provides additional degrees of freedom for radiative transfer, 
the rate of radiative transfer is higher for the tube. 

Figures 5.10 and 5.11 show the velocity and temperature profiles across the channel 
at three different x-locations (x= 2, 5 and 10 cm). It is seen that the flow attains a near 
steady-state value towards the exit of the channel. It is also noted that due to the 
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Fig. 5.6 Radiative flux vs y for H 2 0, NO and OH, P = 1 atm = 3.0. 
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Fig. 5.10 
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boundary layer interaction, the temperature increases along the flow. The radiative heat 
transfer result shows the same trend as Fig. 5.4 

In order to investigate the effect of boundary conditions on the flow, two different 
boundary conditions were selected. In the first case the temperature of the boundary 
was maintained constant at a higher value than the free stream temperature and in the 
second case it was maintained at a lower value. The free stream values were kept 
the same as before. 

Figures 5.12 and 5.13 show the temperature and velocity ptofile for the constant 
temperature wall boundary condition with where T w = 3000 K. It seen that the 
temperature boundary layer is more well developed than is the velocity boundary layer, 
due to energy transfer from the high temperature wall to the fluid. Figure 5.14 shows the 
radiative flux for the above condition. It is seen that near the entrance of the channel, 
due to a high temperature gradient the radiative flux is very high. However, towards 
the exit, since the flow attains a steady-state value, the temperature is evenly distributed, 
resulting in the decrease in radiative flux. 

Figures 5.15 and 5.16 show the velocity and temperature profile when the wall 
temperature is kept lower than the free stream temperature (T w = 1000 K). In this case, 
an opposite trend is seen. Figure 5.17 shows the radiative flux for the above case. The 
result shows the same trend as before. Since there is change in the direction of heat 
flow, the radiative flux changes sign. 

The influence of radiative interactions in chemically reacting supersonic internal 
flows was investigated by considering the physical model shown in Fig. 2.5. The 
specific problem considered is the supersonic flow of premixed hydrogen and air with an 
equivalence ratio of unity in a channel with a compression corner on the lower boundary. 
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Fig. 5.13 Temperature variation across the channel at various x locations, M 
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Fig. 5.14 Variation of q RY with y for different x locations, 
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5.15 Variation of axial 
3.0, T w = 1000 K. 
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Fig. 5.16 








The physical dimensions considered for obtaining results are L=2 cm, x^=l cm, X2=2 
cm, L x =x 1 +X2=3 cm, and a=10°. The inlet conditions, which are representative of the 
scramjet operating conditions, are P=1 atm, T=900 K and M=3.0. The flow is ignited 
by the shock from the compression comer. The flowfield for this problem has been 
investigated by several researchers [25-29], where different models have been used. 
Influence of radiative interactions was investigated in [25] by considering a simple two 
step model. Recently a [28], comparitive study of the flow field was conducted using 
three chemistry models (see Table 2.1). It was found that significant amount of radiation 
participating species are produced by the 35-step chemistry model. Complete discussion 
of the use of the three chemistry models with and without radiation are provided in 
[28], Selected results obtained by using th 18-step chemistry model are presented here 
to demonstrate the influence of radiative interactions. 

The computed results for the 18-step chemistry model are presented in Figs. 
5.18-5.22, both with and without radiative interactions. The results were obtained using 
a 31x31 grid; this was found to be an appropriate grid for the model. The variations in 
temperature, pressure, and species concentration along the x-coordinate are shown for a 
y-location of 0.025 cm from the lower wall. The shock occurs at about x/L x =l/3, and it 
was noted in [28] that the 35-step chemistry model predicts the ignition time accurately. 

The variation in temperature and pressure is shown in Figs. 5.18 and 5.19, 
respectively. The temperature is seen to increase uniformly along the channel (Fig. 
5.18) and there is a significant increase in pressure after the shock (Fig. 5.19). The 
strong fluctuation in pressure after the main shock is attributed to reflection of weak 
shocks from the boundaries. 

Variation in species production are shown in Figs. 5.20-5.22. At moderate 
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temperatures, symmetric molecular species such as O 2 and H 2 do not participate in the 
radiative process. Consequently, under the physical conditions of the problem, mass 
fractions of O 2 and H 2 are not influenced by the radiative interaction (Fig. 5.20). The 
concentration of O 2 and H 2 decreases in the downstream regions as these species react 
to form H 2 O and OH. A slight decrease in the production of both H 2 O and OH is 
noted after the shock due to radiative interactions. In the case of reacting flows without 
radiation, the ignition is seen to take place at x/L x =0.37 indicating a longer ignition 
delay [28]. However, for the reacting and radiating case, the ignition seen to occur at 
x/L x =0.33. This demonstrates that the effect of radiative interaction is to nullify the 
difference in the ignition delay. 

The results of radiative heat transfer by three chemistry models were compared 
in [28], and these are presented here for the sake of completeness. The results for 
the normal radiative flux presented in Fig. 5.23 demonstrate that radiative interactions 
increase rapidly after the shock. The three models are seen to predict the same general 
trend. The results of streamwise radiative flux illustrated in Fig. 5.24 show that the net 
q^x decreases towards the end of the channel. This is due to the cancellation of fluxes in 
the positive and negative x-directions. It is noted that the net radiative transfer is in the 
negative x-direction. The 18-step and 35-step models are seen to predict significantly 
higher q^x than does the 2-step model. This is because radiative heat trransfer is a 
strong function of temperature, pressure, and species concentration, which are higher (in 
the positive x-direction) for the 18-step and 35-step models than for the 2-step model. 
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Fig. 5.22 
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5.23 Variation of normal radiative flux with x for three 
chemistry models. 
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Fig- 5.24 Variation of atreamwise radiative flux with x for three 
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Chapter 6 
CONCLUSIONS 

The two dimensional Navier-Stokes equations have been used to investigate the 
influence of radiative energy transfer on the entrance region flow under supersonic 
flow conditions. In the hypersonic propulsion system, the temperature ranges from 
one to five thousand degrees Kelvin. In this range, various nonsymmetric molecules 
become highly radiative participating. One-dimensional radiative flux was included in 
the energy equation for the solution of nonreacting supersonic flows between parallel 
plates. Specific results have been obtained for different amounts of H 2 O, OH, and NO 
in combination with air. Results demonstrate that the radiative interaction increases with 
increase in pressure, temperature, amount of participating species, plate spacing, and 
Mach number. This can have a significant influence on the overall energy transfer in the 
system. Most energy, however, is transfered by convection in the flow direction. 

The radiative interactions in reacting flows have been investigated by considering 
the supersonic flow of premixed hydrogen and air in a channel with a 10° compression 
corner at the lower boundary. The finite rate chemistry model employed is a nine species 
18-step model. The results show that the effect of radiative interaction is to nullify the 
ignition delay i.e., the reacting and radiating flow predicts a more accurate position of 
the shock. The radiative influence is found to be stronger in the boundary layers, and 
this is seen to change the temperature, pressure, and species concentration in the flow 
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direction. Thus radiative heat transfer can have a significant effect on the entire flow 
field of a hypersonic propulsion system. 

Some earlier results obtained during the course of this study were presented at the 
29 th Aerospace Sciences Meeting in Reno, Nevada, January 7-10, 1991 (AIAA Papers 
91-0373 and 91-0572). These are provided in this report as Appendices F and G. 
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APPENDIX A 


DERIVATION OF THE CONDUCTION HEAT FLUX TERM 

The derivation of the conduction heat flux terms is obtained directly from Ref. 27. 
To simplify Eq. (2— 3a) and Eq. (2-3b) , the Lewis number is assumed to be unity . 
This simplification is carried out in detail for Eq. (2-3a) and the same is applied to 
Eq. (2— 3b) . Using the expression for the thermal diffusivity (a) and Lewis number 
(L e ) , Eq. (2-3a) can be expressed as 

K 

0 = 


L e = 


a 

D 


(A-l) 


n ( T at ^>af iu 

q C y - -pD (^LeCp— + J2 ^ h i 
Defining the binary diffusion coefficient D in terms of the Prandtl and Lewis number 
Eq.( A-l) can be expressed as 


Pr = 


u 

a 


D = — = — - — = — - — 
L e P r Le pP rLe 


qcy — — 


P 

Pr 


- aT 

Cp ay 


m nr 


i=:l 


(A-2) 
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where 


c P - E ^ 

i=l 

The static enthalpy of the mixture is given by the relation 


m 


h = E 


1=1 


h? + / C Pi 


dr] 


fi 


(A-3) 


It should be noted that rj is a dummy variable employed to evaluate the sensible enthalpy 
. Using the Lebnitz formula Eq. (A- 3) is differentiated to obtain 


ah /lO 1 r, , x. \ r ah i P 

fly = El [ i+ J Pl (,) n ) dj + * a^ 


(A-4) 


+ fi / 


aepiW , . fr ,_.aT 

—&T d ” + f s c P, (T ) 


The coefficient on the right hand side is equal to hj and the second and third terms are 
identical to zero , therefore , Eq.(A-4) reduces to 


a , m 

<9h _ ^ 

~ E 

J i=i 


,3 + f.c 
h, ay ‘ "flyj 

Or 


ah _ ^ afj aT 

<9y 11 ay + p ay 

i=l 

Substituting Eq. (A-5) intio Eq. (A-2) q cy is expressed as 

/i ah 


qcy — ~ 


p r a y 


Or 

7 n de 

qcy = 


(A-5) 


(A-6) 
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APPENDIX B 


RADIATIVE FLUX EVALUATION 
FOR THE PSEUDO GRAY MODEL 


Most of the materials in this appendix are also taken from Ref. 27. This appendix 
shows the discretization and evaluation of the radiative flux for the y - direction. The 
same thing can be applied in the streamwise direction. Equations (3-31) and (3-32) 
in the y - direction are written as follows 


d 2 qr(y) 

dy 2 



(y)q r (y) = 3 «(y) 


de(y) 

dy 


(B-l) 


.(y)(i-i) qt (y)l^o-iM |y=0 . 0 


(B-2) 


«(y.|i-i) qr (y)l y=L -iM| y=L =o 


(B-3) 


The above equations is discretized by central differencing . The second derivative of 
q r in the physical domain is discretized as 

d 2 q r . 2 


dy 


2 IJ 


AyjU + tfj) 


qj+i - qj qj - qj-i 
z 3 j A yj A yj 


where 


A yj = yj - yj-i and h = 

yj yj-i 


(B-4) 


Equation (B-4) in discrete forms are written as 
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A 2 yj (1 + /?j) 


^j-i 


1 1 

+ 


9 2 

+ -K 2 


[A yj (l + /3j) V^jAyj Ay jy / 4 Jj 


Hi 


Ay?(i + ^j)/?j 


Hj+i 


= 1.5k: 


e j+l e j , e j e j~l 

^j A yj A yj 


(B- 5 ) 


** 1 6! 2 ) + 3A yi 


qi - 


3 Ayi 


q 2 = 0 


(B- 6 ) 


3A ^j 


7 * 11—1 + 


1 ^ 1 

~ I K i + 


62 2J J 3 Ayj 


Hi 


= 0 


The above can be written in matrix form as 


A 

B 

0 ... 



0 

qr t 

Ri 

C 

A 

B 

0 ... 


...0 


r 2 



0 

c 

A 

B 

qq-x 

= 

0 .. 



0 

C 

A 

q rj 

R j 


where 


(B- 7 ) 


a , 1 l\ 1 

A i,i = *1 77 “ o + 


Aij = 


ei 2/ 3Ay x 

2 / 1 1 \ 9 ? 
A yj C 1 + ^j) v^ Ay J + Ay J/ ^ 

AiJ = Kj (n " 2) + 3Ayj 


; j = 2,j-i 


Bl ,2 = - 


3Ayi 
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B U = 


(i + 0j) h 


; j = 3,j 


c ij = 


Ay?(i + /?j) 


; j = 2,j-2 


dj-1 “ 


3A yj 


Rl = 0 


R; 


15 *j 


e j+l e j , e j e j~l 
^j A yj A ^j 


j = 2, j — 1 


The tridiagonal matrix on the left hand side of Eq. (B-7) is solved efficiently by 
thomas algorithm 

For a circular tube the differential form of radiative flux equation is given by 


d_ 

dr 


1 a . ; 

--r( r qr) 

r dr 


9 2 „ de 

-K q r = 3k— 
4 dy 


«(y) 



2 ) qr|y ^° 3 


1 d r t 
-~r ( r Qr) 
r dr 


ly=0 = 0 




t 2 2 

Differentiating Eq. (B-8) we get 


qr|y=L + 3 


id 

(rqr) 

r dr 


ly=L = 0 


d 2 q r 1 dq r 

dr 2 r dr 





(B-8) 


(B-9) 


(B- 10) 
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Eu. (B-10) along with boundary conditions (B-9) are discretized to give. 
We know by central difference that 


d 2 q r _ 2 

dy 2 ” Ayj(l+/3j) 


qj+i-qj qj-qj-i 

^j A yj A yj 


de 

dr 


1 

2 


e j+i ~ e .i + e j ~ e j-i 


^j A yj 


Ay; 


(B-l 1) 


where 


dq 

dr 


1 

2 


qj+i - qj , qj - 3ti 

/?j A yj Ayj 


Ayj = yj - yj-i 


A = %±LZli 
yj-yj- 1 

Substituting Eqs. (B- 11), into (B-10) we get , 


AyjU+^j) 


qj+i-qj qj-qj-i 
^jAyj Ayj 


- ~K* + 


qj 


(B-12) 


i. 5 KI |%— 2i + 2_2izi) * 


J 1 /3 jAyj 


Ayj 


i / qj+i-qj qj-qj-i 

2r \ /djAyj Ayj 


taking qj+j, qj, qj.j on one side we get, 


+ 


[Ay 2 (1 + /?j)0j Sy^jAyj 


qj+i 


+ 


[Ayj (1 + /5j) \djAyj Ayj 


9 2 1 

4 v" 


qj 


1 

2 


^jAyj Ayj 


qj + 


Ayj(l+/3j) 2yAy 


' i3 i A yj ‘-yj , 

(B- 13) 
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This along with the two boundary condition form a tridiagonal system which can 
be efficiently solved by the Thomas algorithm 
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APPENDIX C 


EXPONENTIAL KERNAL APPROXIMATION 


The Integra — differential governing equation for radiative flux for parallel plate 
geometry, gray non- scattering medium with black bounding surfaces is given by 


q r (y ) = eiuj. 


y 

e 2w +| J F lw-(t) K w ex P 
0 



(t - y) 


dt 




/c^exp 


(t - y) 


dt 


(C-l) 


or it can also be written as 


y 

q r (y) = 2 <tT?F 2 (y) - 2<rTf F 2 (L - y) + 2<r J T 4 (t) k u Fi (y - t) dt 

0 

(C-2) 

L 

-2(7 J T 4 (t) k w F x (t — y) dt 

y 

For Exponential Kemal approxiamation let 


Fi(t) = 


3 


4 


-■2t 

e 2 1 


F 2 (t) = -« w e 


-5t 


(C-3) 
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Substituting Eq. (C-3) into Eq. (C-2) we get 


q r (y) = (7Tf - <rT|e _ f' ,c ‘*'( L - y )+ 


f 

' y 


- L 


2 Kw ) 

/ T 4 (t)e~t^( y_t) dt 

— 

/ T 4 (t)e~3 /c “( t - y )dt 

> 

l 

-0 


J 

Ly J 

j 


(C-4) 


or 


q r (y ) = F (y ) + A | J T 4 (t) e _/? ( y_t )dt - J T 4 (t) e _/? ( t_y )dt J 


where 


F (y ) = aTfe-fr - aT^e~^y~^ 


3 3 

P — t; k u> ; A = /?cr = —K^a 

& 2 


differentiating Eq. (C-4) 


dq r dF 
dy dy 


+ M 


T 4 (y)-0 + 


y 

j T 4 (t), 


-P(y-t) (_ 


(-fldt 


(C-5) 


L 

0 - T 4 (y) + J T 4 (t) e-W-y) (/?) dt 


Simplifying we get ; 
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^7 = + 2At4 (y) + W | - J t 4 (t) e-^y-^dt - J t 4 (t) e -^*-y)dt 

l 0 y 


= ^ + 2at “ (y > - p x 


y l \ 

J T 4 (t)e~^ y_t) dt + J T 4 (t)e- /? ( t - y Mt 1 
0 y J 


(C- 6 ) 


where 


/ 3 A = /? 2 c r = —a 

4 

Differentiating once again we get ; 



d 2 F 

dy 2 


dT 4 

+ 2 A— — j 3 \{ 

dy 


y 

T 4 (y) - 0 + J T 4 (t)e" /? ( y - t )(- / 3 )dt 
0 


+ 


0 - T 4 


L 


(y) + J T 4 (t) (+/?) dt } 


y 


Simplifying Eq. (C- 7 ) results in, 


(C- 7 ) 


d 2 q r 

dy 2 


d 2 F 4 dT 4 


+ ; 3 \ 


$ 


y 

! 


T 4 (t) e' 


-&{y- 


_t )dt -0 J T 4 (t)e~^ t_y) dt J 


(C- 8 ) 
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Eliminating the integrals between Eq. (C-8), one obtains 


d 2 q r 

dy 2 


j2p Jrp4 

d^ + 2A -37 + ' ?2 ^^- F Wl 


or 


d 2 q r 

dy 2 


£ 2 qr(y) = 2A 


dT 4 

dy 


d 2 F 



/? 2 f (y) 


(C-9) 


For the specific case 


Thus 


F(y) = aTfe-fr - aT$e-K L -y'> ; /3 = ^ 


dF 

dy 


-/?<7T|e-^ - op T|e“^ L - y ) 


(C-10) 


d 2 F 

dy 2 


P 2 oT\e~ 0 y - p 2 oT\e~^ L ~^ 


d 2 F 

dy 2 


- p 2 F (y) =0 


(C-ll) 


therefore the governing becomes 


1 d 2 q r (y) 
«w dy 2 


-|qr(y) 


„ O dT 4 
3 J 
dy 


(C- 12) 


Since e = o T 4 


l d 2 q r (y) 9 _ 1 de(y) 

k 2 dy 2 4 qr k u dy 


(C-13) 
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APPENDIX D 


COMPONENTS OF THE JACOBIAN MATRIX 


The material presented in this appendix is taken directly from Ref. 27. In this 
appendix, the species matrix in the left hand side bracket of Eq. (4—9) is evaluated. The 
source term is a function of density, temperature and various species. In evaluating 
the density and temperature dependency is neglected for computational efficiency. The 
components of the Jacobian matrix are as follows : 


Hj = W; = CiMi 
Cj = (pfi/mi) x 10 


_ 3 gm ~ mal 
cm 3 


3W 0j 

5Uq 2 

<9Wq 2 

d U H 2 0 

9Wq 2 

<9U Ha 

gw Qi 

dC 0 H 

5Uo s 


= -KfC Hj 

= o 


= -KfCo, 


= 2K b[ 
= -K fl 


Mq 2 
m oh 
Mh 


OH 


Mo, 


c h 2 


(D-l) 
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-Kf.Co, - K^Cqh 


aw H , 

au H , 

aw H , 

au H ,o 


2K bl 


m h , 

Mh 3 o 


Ch,0 


JWh, 

9Uqh 

Wh,o 

au 0l 

9W H ,0 

9UH a o 


Mh Mu 

= 2K b ,^C OH -2K fj — it-CoHC Hl 

Mqh m OH 

= 0 


= -4K b2 C H3 o 


^w Hi o 

dy H , 

d^H,,0 

5Uqh 


= 2K fj 
= 4K f2 


M H 2 0 

Mh 2 

m h 2 o 

Mh 2 


c 2 

^OH 

C OH C H 2 


^OH 

5Uo 2 

dW OH 

3UhjO 


2Kf ^2flc H 


4Kk 


Mq 2 
m oh 


M 


h 2 o 


h 2 o 


^OH 

<9Uh 2 

<3Wqh 

5U OH 


2K 
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fi 


OH, 


M 


h 2 


o 2 


— 2K 


m oh 

f! M Hl 


c 


2 

OH 


-4K b2 CpH - 4Kp 2 CpjjC H2 


(D-2) 
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APPENDIX E 


PROGRAM FOR THE CALCULATION OF 
PLANCK MEAN ABSORPTION COEFFICIENT 


c ********************************************************** 
c Program To Calculate Kappa ( Kp ) For different gases 
c Provision is made for calculating Kp for H20, OH, NO 
c Various combination of these gases can be taken and Kp 
c can be determined 

c********************************************************** 

Program Kappa 

parameter (ix=31, iy=31, iq=12, is=9, ip=10, im=8, ir=18) 
common/ rad/ omega (7) , ski (ix, iy) , sk2 (ix, iy) , phi 7 (ix, iy) , 
&eo(ix,iy,7) ,co2 (ix, iy, 7) , rkpt (ix, iy) , 

&ao (ix, iy, 7) , uc (ix, iy, 5) , rkp (ix, iy, 7) , cc (52, 52) , 

&aa (52 , 52 ) ,bb(52,52) ,b2 (ix, iy, 5) , et (ix, iy ) ,btn(ix,iy) , 
&psp (ix, iy, is) , rsp (5) ,dd(52,52) , qrx (ix, iy) , qry (ix, iy) , 
&rkpt 1 (ix, iy) , rkpt 2 (ix, iy) , phiv (ix, iy, is) , fn (ix, iy, is) 
dimension estop (ix,iy) 

dimension irror (100, 2) , an (31, 31 ) , bn (31, 31) , cn ( 3 1 , 31) 
dimension dn ( 31 , 31 ) , r (31 , 31 ) , pn (31,31) , rhon (31, 31 ) 
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integer btl, bt2, bt3, bt4 

open (unit=26, f ile= ' trash' ) 

do 1 izz = l , 31 

r (izz, 1) =8314 .34 

btn (izz, 1) =izz*150 

fn (izz, 1, 1) =0 . 0 

fn ( izz, 1 , 2 ) =0 . 0 

fn (izz, 1, 3) =0 . 99 

fn (izz, 1, 4) =0 . 00 

fn (izz, 1, 5) =0 . 0 

fn (izz, 1 , 6) =0 . 0 

fn(izz,l,7)=0.00 

fn (izz, 1, 8) =0 . 00 

fn (izz, l,ncs) =0 . 01 

pn (izz, 1) =101325 . 

rhon (izz, 1) =pn(izz, l)/(r(izz,l) *btn(izz, 1) ) 
1 continue 
nnx=ix 
nn=iy 
nxpl=ix-l 
nypl=iy-l 
ix j=ix 
iy j=iy 
izz=31 

nxpl = nnx+1 
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nypl = nny+1 


c 

c initialize omega (wave no. for planck' s formula) 

c at band center 

c 


c 

c 

c 


omega (1) = 
omega (2 ) = 
omega (3) = 
omega (4) = 
omega (5) = 
omega ( 6 ) = 
omega (7) = 
soh = 110. 


500. 

1600. 

3750. 

5350 . 

7250. 

3570. 

1876 

* 10000 . /101325 . 


initializing for band absorption 


eh = 

6 . 625e-27 

c = 

2 . 998e+10 

urk 

= 1 . 380e-16 

ggl 

= 3652. 

g2 = 

1595. 

g3 = 

3756. 

cl = 

3 .1415926 

c2 = 

eh * c/urk 

vl = 

1 . 
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v2 = 0. 
v3 = 1 . 
epsi = 1. 
depsi = 1 . /epsi 
do 88890 izz=l,ixj 

skl(izz,l) = (btn ( izz, 1 ) /300 . ) **0 . 5 
sk2(izz,l) = (300 . / (btn (izz, 1) ) ) **1 . 5 
b2(izz,l,l) = . 073/skl (izz, 1) 
b2(izz,l,2) = . 130/skl (izz, 1) 
b2(izz,l,3) = . 145/skl (izz, 1) 
b2(iz'z,l,4) = . 118/skl (izz, 1) 
b2(izz,l,5) - . 201/skl (izz, 1) 

88890 continue 

do 88960 izz=l, ixj 

phi7(izz,l) = (btn ( izz , 1 ) /100 . ) * * (-0 . 5) 
phi7(izz,l) = -17.6 * phi7(izz,l) 
phi7(izz,l) — exp (phi7 ( iz z , 1 ) ) 

88960 continue 

do 88990 izz=l,ixj 

phiv ( izz , 1 , 1 ) = (-eh*c*ggl) / (btn (izz, 1) *urk) 

phiv ( izz , 1 , 1 ) = exp (phiv (izz, 1, 1) ) 
phiv ( iz z , 1 , 1 ) = 1 . -phiv (izz , 1 , 1 ) 

88990 continue 

do 89020 izz=l, ixj 

phi v ( izz , 1 , 2 ) = (-eh*c*g2) / (btn (izz, 1) *urk) 
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phi v (iz z , 1 , 2 ) = exp (phiv (izz, 1, 2) ) 
phiv (izz, 1, 2) = 1 . -phiv ( izz, 1 , 2 ) 

89020 continue 

do 89050 izz=l,ixj 

phiv ( izz , 1 , 3 ) = (-eh*c*g3) / (btn (izz, 1) *urk) 

phiv ( izz , 1 , 3 ) = exp (phiv (izz, 1, 3) ) 
phiv ( izz , 1 , 3 ) = 1 . -phiv ( izz, 1 , 3) 

89050 continue 
c 

c calculate denominator of phi-101 

c 

do 89080 izz=l,ixj 

phiv ( izz , 1 , 1 ) = (phiv (izz, 1, 1) * phiv (izz , 1 , 2 ) 

& * phiv ( izz , 1 , 3) ) 

89080 continue 
c 

c calculate numerator of phi-101 

c 

do 89090 izz=l,ixj 

phi v ( izz , 1 , 2 ) = (-eh*c* (vl*ggl+v2*g2+v3*g3) ) / 

& (btn ( izz , 1 ) *urk) 

phiv ( izz , 1 , 2 ) = exp (phiv ( izz, 1 , 2 ) ) 
phiv ( izz , 1 , 2 ) = 1 . -phiv ( izz , 1 , 2 ) 

89090 continue 
c 
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c 


calculate phi-101 


do 89120 izz=l,ixj 

phiv { izz , 1 , 3 ) = phiv ( izz , 1 , 2) /phiv ( izz , 1 , 1 ) 

89120 continue 

do 2 k=l , 7 
do 89130 izz=l,ixj 

eo(izz,l,k) = c2 * omega(k) /btn(izz,l) 
eo(izz,l/k) = exp (eo (izz , 1 , k) ) 

eo ( izz, 1 , k) = cl * {omega (k) **3 .)/ (eo (izz , 1 , k) -1 . ) 
eo(izz,l,k) = eo(izz,l,k) * 1.0e-05 
89130 continue 
2 continue 

do 89170 izz=l,ixj 

co2 (izz, 1,1) = 771. * sk2 ( izz , 1 ) * phi7(izz,l) 

co2 (izz, 1,2) = 3.35 * sk2(izz,l) 
co2(izz,l,3) = 1.52 * sk2 (izz, 1) 

co2(izz,l,4) = 0.276 * sk2(izz,l) * phiv (izz, 1, 3) 
co2(izz,l, 5) = .230 * sk2(izz,l) * phiv (izz, 1, 3) 
89170 continue 
c 

c unit change 

c 

do 9 k=l,5 
do 89220 izz=l, ixj 
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co2 ( izz , 1 , k) = (100 . /101325 . ) * co2(izz,l,k) 
co2(izz,l,k) = co2(izz,l,k) 

89220 continue 
9 continue 
izz = izz 

do 89230 izz=l,ixj 

if ( . not . ( f n (izz + 1 , 2 , 3) . le . 0 . ) ) go to 89230 
f n ( izz+1 , 2, 3) = 0.0001 
89230 continue 
izz = izz 


do 89240 izz=l,ixj 


psp (izz, 1, 3) 
& *btn (izz , 1 ) 
psp (izz, 1,3) 
psp (izz, 1,4) 
& *btn (izz , 1 ) 
psp (izz, 1, 5) 
& *btn (izz, 1) 
89240 continue 


rhon (izz, 1) *fn(izz,l,3) *r(izz, 1) 


1 . 

rhon (izz, 1) *fn (izz, 1, 4) *r (izz, 1) 
rhon (izz, 1) *fn (izz, 1, 8) *r (izz, 1) 


c 

c calculate ao 

c 


do 89260 izz=l,ixj 
ao(izz,l,l) = 49.4 * skl(izz,l) 
ao (izz, 1,2) = 90.1 * skl(izz,l) 
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ao(izz,l,3) = 112.6 * skl(izz,l) 
ao(izz,l,4) = 79.7 * skl(izz,l) 
ao(izz,l,5) = 79.7 * skl(izz,l) 

89260 continue 
c 

c change unit 

c 

do 4 k=l,5 
do 89310 izz=l,ixj 
ao(izz,l,k) = ao(izz,l,k) * 100. 
ao(izz,l,k) = ao(izz,l,k) 

89310 continue 
4 continue 
c 

c Calculating Kp for various gases and their combination 

c 

do 89320 izz=l, ixj 

et(izz,l) = 5 . 668e-05*btn (izz, 1) **4 . 
rkpt (izz , 1) = 0 . 

89320 continue 

do 3 k=l , 5 
do 89340 izz=l,ixj 

co2 (izz, 1, 6)=110.*300./ (10*btn (izz, 1) ) 
ao (izz, 1 , 6) =1 . 

co2 (izz, 1,7)=132.*300./ (10 . *btn (izz, 1) ) 
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ao (izz, 1, 7) =1 . 

rkp(izz,l,k) = eo (izz , 1 , k) *ao (izz, 1 , k) *co2 ( izz , 1 , k) 
rkpt (izz, 1) = rkpt(izz,l) + rkp(izz,l,k) 
rkpt 1 (izz, 1 ) =eo (izz, 1, 6) *ao (izz, 1, 6) *co2 (izz, 1, 6) 
rkpt2 (izz, 1) =eo (izz, 1,7) *ao (izz, 1,7) *co2 (izz, 1,7) 
89340 continue 
3 continue 

do 89360 izz=l, ix j 

rkpt (izz, 1) = rkpt (izz, 1) *psp (izz, 1, 3) /et (izz, 1) 
rkptl(izz,l) = rkpt 1 (izz, 1) *psp (izz, 1 , 4 ) /et (izz, 1 ) 
rkpt2 (izz, 1) = rkpt 2 (iz'z , 1 ) *psp (izz, 1 , 5 ) /et ( izz , 1 ) 
rkpt ( izz , 1 ) =rkpt (izz, 1 ) + rkpt 2 (izz , 1) 
write (26, 20) btn (izz, 1) , rkpt (izz, 1) 

20 format (2x, f 10 . 4, 2x, e20 . 6) 

89360 continue 
stop 
end 
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RADIATIVE INTERACTIONS IN A HYDROGEN-FUELED SUPERSONIC COMBUSTOR 

R. Chandrasekhar* and S. N. Tiwari* 

Old Dominion University. Norfolk, VA 23529-0247 
and 

J. P. Drummond 1 
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Abstract 

The two-dimensionaJ, elliptic Navicr-Stokcs equa- 
tions are used to investigate supersonic flows with 
finite-rate chemistry and radiation, for hydrogen-air 
systems. The chemistry source term in the species 
equation is treated implicitly to alleviate the stiffness 
associated with fast reactions. The explicit, unsplit 
MacCormack finite-difference scheme is used to ad- 
vance the governing equations in time, until conver- 
gence is achieved. The specific problem considered is 
the premixed flow in a channel with a ten-degree com- 
pression ramp. Three different chemistry models are 
used, accounting for increasing number of reactions 
and participating species. Two chemistry models as- 
sume nitrogen as inert, while the third model accounts 
for nitrogen reactions and NO, formation. The tan- 
gent slab approximation is used in the radiative flux 
formulation. A pseudo-gray model is used to repre- 
sent the absorption-emission characteristics of the par- 
ticipating species. Results obtained for specific condi- 
tions indicate that the radiative interactions vary sub- 
stantially, depending on reactions involving HQ* and 
NO species, and that this can have a significant influ- 
ence on the flowficld. 

Nomenclature 

A band absorplance , m -i 
,-l„ band width parameter, m~* 

C } concentration of the j* h secies, ky-mole/m 3 
(' v constant pressure specific heat , J/kg - !\ 

('u correlation parameter , (Af/m 2 )~ l ii*“ l 
h' total internal energy 
Planck's function 
fj mass fraction of j 1 * 1 species 
II total enthalpy, J /leg 
h sialic enthalpy, J/kg 
k thermal conductivity 
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kb backward rate constant 
kj forward rate constant 
P pressure , N/m 7 
Pj partial pressure of j th species 
<jn total radiative flux’ 

R gas constant 

S integrated band intensity, (jV/m 2 )“ l iu~- 

T lemjierature , l\ 

u, v velocity in x - and y - directions , m/s 
Wj production rate of j th species, kg/nP - s 
x » y physical coordinates 
7 ratio of specific heats 
Kf, Planck mean absorjtt ion enc f f ictt ut 
A second coefficient of viscosity , wan h nyih 
P dynamic Viscosity , kg/m - s 
£> V computational coordinates 
p density 

a Stefan — Boltzmann constant 
t shear stress 
<j> equivalence ratio 
w ware number, m ~ 1 

Introduction 

In recent years there has been a renewed 
interest in the development of a hypersonic transat- 
mospheric aerospace vehicle capable of flying at sub- 
orbital speeds. A hydrogen-fueled supersonic combus- 
tion ramjet (scramjct) engine is a strong candidate for 
profiling such a vehicle. For a better understanding 
ot the complex flowficld in different regions of the en- 
gine, both experimental and computational techniques 
arc employed. Several computer programs have been 
developed 1 1—4] and applied to gain more insight into 
the problem involving the flow in the various sections 
of the scramjct module. 

The combustion of hydrogen and air in the scram - 
jet combustor results in absorbing-emitting gases such 
as water vapor and hydroxyl radicals. Existence of 
such gases makes llie study of radiation heat transfer an 
important issue. There are several models available in 
the literature to represent die absorption-emission char- 
acteristics of molecular gases 15-10). One- and two- 
dimensional radiative heat transfer equations for vari- 
ous flow and combustion related problems arc available 
1 1 1- 19(. In earlier studies 1 16,18,19), both pseudo-gray 
and nongray gas mtxlcls were employed to evaluate ra- 
diative heal transfer for chemically reacting supersonic 
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flow. Results of both models were compared and the 
pseudo-gray model was found to be computationally 
more efficient 

Considerable work has been carried out in the past 
decade to model the chemical kinetic mechanism of the 
hydrogen-air system. A most complete model would 
involve some 60 reaction paths [20], rendering numeri- 
cal solution very difficult, if not impossible. A two-step 
chemistry model, has been used for computing super- 
sonic combustion [4, 16, 18, 19]. This model has only 
four species and two reaction paths, and is useful for 
preliminary studies. However, there are several limita- 
tions to this model, such as ignition-phase inaccuracy 
(i.e. a much shorter ignition delay) and also, overpre- 
diction of flame temperature and longer reaction times. 
Recent improvements in this area include a 8-species, 
14-reaction model [21] and a 9-species, 18-reaction 
model [2, 22]. While none of these aforementioned 
models account for nitrogen reactions (assuming ni- 
trogen as inert), recent developments in this area in- 
clude a 15-species, 35-rcaction model which accounts 
lor NO, formation and other nitrogen reactions in the 
hydrogen-air system [22]. 

The objectives of the present study are to extend 
the radiative heat transfer formulation used with the 
global two-step chemistry model [18, 19], to the more 
complete models namely the 9-species, 18-reaction 
model as well as the 15-species, 35-reaction model. 
The effect of radiative heat transfer in both transverse 
and streamwise directions is investigated. The finite- 
difference method using the explicit, unsplit MacCor- 
mack scheme [23] is used to solve the governing equa- 
tions. 

The flowficld in the combustor is represented by 
the Navier-Stokcs equations and by the appropriate 
species continuity equations [2, 3]. Incorporation of 
the fmite-rale chemistry models into the fluid dynamic 
equations can create a set of stiff differential equations. 
Stillness is due to a disparity in the time scales of the 
governing equations. In the Lime accurate solution, af- 
ter the fast transients have decayed and the solutions 
are changing slowly, taking a larger time step is more 
efficient. But explicit methods still require small time 
steps to maintain stability. One way around this prob- 
lem is to use a fully implicit method. However, this 
requires the inversion of a block multi-diagonal system 
of algebraic equations, which is also computationally 
expensive. The use of a semi-implicit technique, sug- 
gested by several investigators [24-26], provides an 
alternative to the above problems. This method treats 
the source term (which is the cause of the stiffness) 
implicitly, and solves the remaining terms explicitly. 

Basic Governing Equations 

The physical model for analyzing the flowficld in 
a supersonic combustor is described by die Navier- 


Siokes and species continuity equations. For two- 
dimensional flows, these equations are expressed in 
physical coordinates as, 


dU 9F 0G 

aT + a7 + ^ + // = 0 


(i) 


where vectors U, F, G and H are written as. 
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The viscous stress tensors in the F and G terms are 
given as, 


_ . / flu dv\ On 

A Ur + 0y) ~ 2l ‘{h 

(Du Dv\ 

r - = Ar* + Ty) 

/flu dv\ dv 

~{dl + d^)~ 2, % 


(2a) 

(2b) 

(2c) 


where A = -*^/i . The quantities q cx and q cy in the 
F and G terms arc the components of the conduction 
heat flux and are expressed as 


3T 

qct ~ ~ k HI 

0T 

7cy ' Dy 


(3) 


The molecular viscosity p is evaluated from the 
Sutherland's formula. The total internal energy E in 
Eq. (2) is given by 


_ l‘ u* + t' 2 
— -j- 
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Specific relations are needed for the chemistry and ra- 
diative flux terms. These are discussed in the following 
sections. 


The gas constant for the mixture is evaluated by 
a mass-weighted summation over all species as 


Chemistry and Thermodynamic Model 
Chemical reaction rate expressions are usually de- 
termined by summing the contributions from each rel- 
evant reaction path to obtain the total rate of change 
of each species. Each path is governed by a law of 
mass action expression in which the rale constants can 
be determined from a temperature dependent Arrhenius 
expression. The reaction mechanism is expressed in a 
general form as 

n* k fi< 

f - *'*»."*■ (5) 

> = 1 S > = 1 


where ns = number of species and nr ■ number of 
rcaclions. The chemistry source terms Wj in Eq. (1) 
are obtained, on a mass basis, by multiplying the molar 
changes and corresponding molecular weight as 


nr 

M i C i = 

» = 1 

r a 

- *». n c " r 


n Cm" j 


fft — 1 


»n = l 

• j - 1. 


( 6 ) 


The reaction rate constants kr , and k(, , appearing 
in Eqs. (5) and (6) are determined from an Arrhenius 
rate expression as 


R = 02 ) 

i=t 

P = PRT (13) 


Radiation Transfer Model 

Evaluation of the energy equation presented in Eq. 
(1) requires an appropriate expression for the radiative 
flux term, q R . Therefore, a suitable radiative transport 
model is needed. Various models are available in the 
literature to represent the absorption-emission charac- 
teristics of the molecular species [10]. The equations of 
radiative transport are expressed generally in integro- 
differential forms. The integration involves both the 
frequency spectrum and physical coordinates. In many 
realistic three-dimensional physical problems, the com- 
plexity of the radiative transport equations can be re- 
duced by introduction of the tangent-slab approxima- 
tion. This approximation treats the gas layer as a one- 
dimensional slab in evaluation of the radiative flux (Fig 
1 ). 

Detailed derivations of radiative flux equations for 
gray as well as nongray radiation have been carried out 
previously [15, 19). For a mulliband gaseous system, 
the nongray radiative flux in the normal direction is 
expressed as 


k 


J, 



(7) 


where , 



The coefficients A , N and E appearing in Eq. (7) 
are given in Table 1. The term An in Eq. (9) denotes 
the difference in the number of moles of reactants and 
products. 

The Gibbs energy term in Eq. (9) is calculated as 


9«(y) = «i - e 2 + 


P-'!u 


; 3ll 0 , 

Ai 2~L (y ~ + 


L 

f f<ku,.(;)l , 3 u 0l , 

/ [— 


(14) 
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The information on the band absorptancc A, and other 
quantities is available in the cited references. 

For a gray medium, the spectral absorption coef- 
ficient is independent of the wave number, and 
an expression for the radiative flux is obtained as [5, 
16. 19] 


H • J - 1. nr (10) 
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It is computationally mote efficient to use Eq. (IS) 
in the general energy equation than Eq. (14). This is 
because by differentiating Eq. (IS) twice (using the 
Leibnitz formula) ihe integrals are eliminated and the 
following inhomogenous ordinary differential equation 
is obtained : 


i <* 3 ?«(y) 
** dy a 



3 dejy) 
k dy 


(16) 


The solution of Eq. (16) requires two boundary con- 
ditions which are given for non-black diffuse surfaces 
as [5} 


(.7 ^ y)J v= t 



d<iR 

<*y . 


y = L 


=S 0 
= 0 


(17a) 

(17b) 


For black surfaces t\ « » 1 and Eqs. (17) reduce 

10 simpler forms. 

An appropriate model for a gray gas absorption 
coefficient is required in Eqs. (15) — (17). This is 
represented by the Planck mean absorption coefficient, 
which is expressed for a multi-band system as [5, 19] 


K = Kf> 


Pj 

*T*(y) 


^e Ui (T)S,(7') 
1 = 1 


(18) 


It should be noted that k? is a function of the tem- 
perature and the partial pressures Pj of the species. 


Method of Solution 


The governing equations are transformed from the 
physical domain (x , y) to a computational domain (( , 
?/), using an algebraic grid generation technique similar 
to the one used by Smith and Weigel (27). In the 
computational domain, Eq. (1) is expressed as 


dO_ 5F dG 
di + di + du 


+ fi = o 


(19) 


where 

0 = UJ , F = Fy n - Gx n 
6 = Gx ( - Fy { , H = II J (20) 

J = x^,, - y^x n 


Once the temporal discretization has been performed, 
live resulting system is spatially differenced using 
live explicit, unsplit MacCormack predictor-corrector 
scheme [23j. This results in a spatially and temporally 
discrete, simultaneous system of equations at each grid 
point [25, 26]. Each simultaneous system is solved, 
subject to initial and boundary conditions, by using 


the Householder technique [28, 29]. At the super- 
sonic inflow boundary, all flow quantities are speci- 
fied as frees tr earn conditions. At the supersonic out- 
flow boundary, non-refleclivc boundary conditions arc 
used, i.e. all flow quantities are extrapolated from in- 
terior grid points. The upper and lower boundaries are 
treated as solid walls. This implies a non-slip boundary 
condition (t.e. zero velocities). The wail temperature 
and pressure are extrapolated from interior grid points. 
Initial conditions are obtained by specifying freestream 
conditions throughout the flowfieid. The resulting set 
of equations is marched in time, until convergence is 
achieved. The details of the radiative flux formulation 
and method of solution are available in [19]. 

Results and Discussion 

Based on the theory and computational procedures 
described previously, an algorithm has been developed 
to solve the two-dimensional Navicr-Stokes equations 
for chemically reacting and radiating supersonic flows. 
The extent of radiative heat transfer in supersonic flows 
undergoing hydrogen-air chemical reactions, has been 
investigated using three chemical kinetics models, ac- 
counting for increasing number of reactions and par- 
ticipating species. For the temperature range consid- 
ered in this study, the important radiating species are 
OH and HjO. The gray gas formulations arc based 
on the Planck mean absorption coefficient which ac- 
counts for the detailed information on different molec- 
ular bands. The radiative fluxes have been computed 
using this ‘pseudo-gray’ formulation. The justification 
for using this model is provided in [19]. 

The specific problem considered is the supersonic 
flow of premixed hydrogen and air (stoichiometric 
equivalence ratio <f> » i.O ) in a channel with a com- 
pression comer on the lower boundary (Fig. 2). The 
physical dimensions considered for obtaining results 
arc L = 2 cm., X t * 1cm., X 2 = 2cm., L* » Xi + 
Xz = 3 cm. , and <* = 10 degrees . The flow is ignited 
by the shock from the compression comer. The inlet 
conditions which are representative of scramjct oper- 
ating conditions, are Poo = 1.0 atm., T^ = 900 K 
and Mo, = 4.0. This same flow has been computed 
by several CFD research groups [4, 18, 19, 21 1 as a 
benchmark case. 

Figures 3-6 show the computed results using a 31 
x 3 1 grid , for temperature and pressure as well as H 2 O 
and OH species mass fractions, varying along x at the 
location y = 0.02 cm from the lower wall (boundary 
layer region). Figures 3 and 4 show the temperature 
and pressure profiles predicted by the three chemistry 
models. The temperatures in the boundary layer show 
a gradual increase(Fig. 3). The pressure profiles are 
plotted at y * 0.13 cm. (inviscid region) and show a 
sharp increase due to the shock (Fig. 4). The ignition- 
phase inaccuracies of the three chemistry models can 
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be seen in Figs. 5 and 6. The shock is occurring after 
x / Lx = 0.3. However, the 2-step model predicts 
ignition before the shock (shorter ignition delay) due 
to the high temperature in the boundary layer. On the 
other hand, the 18-step model predicts a longer ignition 
delay, at x / L x - 0.37 (Fig. 5). The 35-step model’s 
prediction of ignition delay appears to be an average of 
the other two models. Although the three models do not 
differ much in prediction of temperature and pressure 
profiles, they do differ significantly in predictions of 
species productions (Figs. 5, 6). 

In order to resolve this discrepancy, a grid sen- 
sitivity study was carried out to examine whether the 
grid size affects the flow predictions. The results of 
three grid distributions 11 x 31 , 31 x 31 and 61 x 
31 are shown in Fig. 7, and it appears that the 31 x 
31 grid is sufficient for the present study. 

The reason for the varying predictions of species 
production by the three models was further examined 
and the results are shown in Figs. 8 and 9. Figure 8 
shows that the Reaction No. 8 in Table 1 is critical 
in determining the extent of chemical heat release and 
H 2 0 production. Reaction No. 8 deals with produc- 
tion of H() 2 radical. This reaction is absent from the 
2-slep model, while it is common to both the 18-step 
and 35-step models. Figure 8 shows that the 35-step 
model experiences nearly a 30% drop in temperature 
at the channel exit, when the rate of Reaction No. 8 
is reduced by a factor of 1000 (effectively culling off 
the production of the H0 2 radical). In contrast, the 
18-step model shows a 15% drop in temperature, when 
subjected to the same reduction in rate of Reaction No. 
8. This shows that the Reaction No. 8 controls dte 
overall H 2 0 production occurring in Reaction Nos. 
9-18 (Table 1). Due to the high temperatures ( " 3000 
K) in the flowfield, there is a pool of highly reactive 
free radicals like H , O, etc. The H0 2 radical is 
converted to the very reactive OH radical, by the free 
radicals (Reaction Nos. 11 and 12). This establishes 
die H0 2 radical as a very important species in pro- 
moting tlame propagation in hydrogen-air flames. A 
similar study has been carried out in (30). Since the 
2-siep model docs not have the H0 2 radical, it pre- 
dicts lesser amounts of OH and H 2 0. 

It was necessary to determine the reason for the 
higher sensitivity of the 35-step model to the H0 2 rad- 
ical, as compared to the 18-step model. Figure 9 shows 
that the Reaction Nos. 2 1 and 23 in Table 1 are critical 
iii determining the extent of chemical heat release and 
ll 2 0 production. Reaction Nos. 21 and 23 deal with 
production of the NO radical. These reactions are 
absent from the 2-siep and 18— step models, whereas 
they play an important role in the 35-step model. Fig- 
ure 9 shows that the 35-step model undergoes a 30% 
reduction in temperature, when the rates of Reaction 
Nos. 21 and 23 arc reduced by a facmr of I (XX) (ef- 


fectively cutting off the production of the NO radical). 
This is nearly the same reduction caused by reducing 
the rate of Reaction No. 8 by a factor of 1000. Due 
to the high temperatures in the flowfield, the usually 
inert nitrogen dissociates into the highly reactive N 
free radical. This free radical N is then oxidized in 
Reaction Nos. 21 and 23, thereby producing the NO 
radical. This NO radical converts the HO J radical 
into the highly reactive OH radical, through Reaction 
No. 29. This confirms that the NO radical is a very 
important species for flame propagation in a hydrogen- 
fueled supersonic combustor. Since the 35-step model 
has the NO radical, it predicts higher amounts of OH 
and H 2 0 than the 18-step model. 

Based on the above understanding of the chemi- 
cal kinetics of supersonic hydrogen-air flames, the ra- 
diative interactions were examined. Figure 10 shows 
the profiles of the normalized sireamwise radiative flux 
qiu predicted by the three chemistry models, along 
the location y * 0.02 cm. from the lower wall. The 
qiu flux reduces towards the end of the channel due 
to cancellation of fluxes in positive and negative direc- 
tions. It is seen from Fig. 10 that the 18-stcp and 
35-slcp models predict significantly higher amounts 
(50% more and 100% more , respectively) of q*. than 
the 2-stcp model. This is because radiative heat tran- 
fer is a strong function of temperature, pressure and 
species concentrations. So the larger values of radiative 
fluxes are caused by higher amounts of H 2 0 concen- 
trations, which in turn, depend on reactions involving 
H0 2 and NO species. 

Figure 1 1 shows the variations of the normal ra- 
diative flux qn y along x , at the location y = 0.02 cm. 
These do not appear to vary significantly between the 
three chemistry models. However, in all three cases, 
the qit y value increases rapidly after the shock. 

Figures 12-15 show the computed results for re- 
acting flows with and without radiation, for the three 
chemistry models. It is seen that the 2-stcp model 
shows only slight effect of radiative interaction, as 
compared to the 18-stcp and 35-step models. The 
18-step and 35-slcp models predict lower tempera- 
ture and lower H 2 0 and OH concentrations after the 
shock. This is because of the q K , flux, which reduces 
the total energy. Comparison of results in Figs. 12-15 
shows that the 35-slcp model exhibits stronger clfcct 
of radiative interactions, than the other two models. 

For reacting flows without radiation, it was seen 
from Figs. 5 and 6 that the 1 8-stcp model hail a longer 
ignition delay (ignition at x / L x = 0.37), while the 
35-step model had a shorter ignition delay (ignition at 
x / L x = 0.27). Another effect of radiative interactions, 
seen in Fig. 14, is to nullify this difference in predic- 
tions of ignition delay. For both 18-stcp and 35-step 
models, widi radiation, die ignition is seen to occur at 
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the same point, x / Lx ■ 0.33. No such effect is seen 
on the ignition characteristics of the 2-step model 

Conclusions 

The two-dimensional, spatially elliptic Navier- 
Stokes equations have been used to obtain solutions 
for supersonic flows undergoing finite-rate chemical 
reactions along with radiative interactions. The spe- 
cific problem considered is of the premixed flow in a 
channel with a ten-degree compression ramp. The in- 
let conditions used in the present study correspond to 
typical flow conditions of a scram jet engine. Three 
different chemistry models were used for parametric 
studies, accounting for increasing number of reactions 
and participating species. It is seen that the radiative in- 
teractions vary significantly, depending particularly on 
chemical reactions involving HO 2 and NO species. 
These reactions have a substantial effect on the flow- 
field, with regard to H 2 0 concentration , temperature 
and pressure. Also, it is observed that the difference 
in the ignition delays of two chemistry models invov- 
ing H0 2 reactions is nullified as a result of radiative 
interaction. The results also show that the strcumwisc 
radiative flux reduces the temperature and concentra- 
tion of species. This effect is a strong function of the 
amount of H 2 0 species concentration. 
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Abstract 

Analyses and numerical procedures are presented to investigate the radiative interactions 
of absorbing-emitting species in chemically reacting supersonic flow in various ducts. Specific 
attention is directed in investigating the radiative contributions of H 2 0, OH, and NO under 
realistic physical and flow conditions. The radiative interactions in reacting flows are investigated 
by considering the supersonic flow of premixed hydrogen and air in a channel with a compression 
corner at the lower boundary. The results indicate that radiation can have significant influence 
on the flowfield and species production depending on the chemistry model employed. 

Nomenclature 


A band absorptance, m' 1 

A () band width parameter, m' 1 

Cj concentration of the j th species, kg-mole/m 3 

e. . Planck’s function 

E total internal energy 

fj mass fraction of j** 1 species 

h static enthalpy, J/kg 

k thermal conductivity, J/m-s-k 
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kb 

kf 

P 

Pj 

qr 

qRa; 

R 

S 

T 

u,v 


x,y 

K(jj 

Kp 

A 

P 


^r, 


P 

a 

T 


ijJ 


backward rate constant 
forward rate constant 
pressure, N/m 2 
partial pressure of j* species 
total radiative heat flux, J/m 2 -s 
spectral radiative heat flux, J/m 3 -s 
gas constant, J/kg-K 

integrated band intensity, (N/m 2 )'^' 2 
temperature, K 

velocity in x- and y-direction, m/s 
production rate of j* species, kg/m 3 -s 
physical coordinates 


spectral absorption coefficient, m ' 1 

Planck mean absorption coefficient 

✓ 

second coefficient of viscosity* wave length, mj 

A 

dynamic viscosity (laminar flow), kg/m-s 

computational coordinates 

density 


Stefan-Boltzmann constant, erg/S-cm 2 -K 4 
shear stress 
equivalence ratio 
wave number, m 1 


Introduction 

There is a renewed interest in investigating various aspects of radiative energy transfer in 
participating mediums. Radiative interactions become important in many engineering problems 
involving high temperature gases. Recent interest lies in the areas of design of high pressure 
combustion chambers and high enthalpy nozzles, entry and reentry phenomena, hypersonic 
propulsion, and defence oriented research. 



Basic formulations on radiative energy transfer in participating mediums are available in 
standard references [1-5]. The review articles presented in [6-15] are useful in understanding 
the radiative properties of participating species and the nature of nongray radiation. The validity 
of radiative transfer analyses depends upon the accuracy with which absorption-emission and 
scattering characteristics of participating species are modeled. There are several models available 
to represent the absorption-emission characteristics of molecular species and these are reviewed 
in [12, 13]. These models have been used to investigate radiative interactions in several duct 
flows [16-29], 

The purpose of this study is to investigate the effect of radiative heat transfer in chemically 
reacting supersonic flow in various ducts under different physical and flow conditions. This pro- 
vides essential information for investigating the effect of radiative interactions in the combustor 
of a supersonic combustion ramjet (scramjet) engine. This hydrogen-fueled engine is proposed 
for propelling transatmospheric hypersonic vehicles. Several basic codes have been developed 
to compute the flowfield in a scramjet engine [21-24]. The combustion of hydrogen and air 
results in absorbing-emitting gases such as H 2 O, OH, and NO. Specific attention, therefore, is 
directed in investigating the radiative contributions of these gases under realistic conditions. In 
essence, the present effort is a continuating of the earlier work conducted in this area of re- 
search [25-27]. Extensive literature survey is provided in the cited references. A comparison of 
different chemistry models used in investigating radiative interactions is presented in [28]. 

Basic Governing Equations 

The physical problem considered to investigate the effect of radiative interactions in super- 
sonic flow are two-dimensional laminar flow between two parallel plates (Fig. la) and within a 
circular tube (Fig. lb). Another geometry is considered to study the effect of shocks and chem- 
ical reactions on the radiative heat transfer and this consists of a channel with a compression- 
expansion ramp (Fig. lc). The governing equations and boundary conditions are provided here 
for all physical problems considered in this study. 
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The physical problem considered for basic understanding of radiative interaction in super- 
sonic flows is two-dimensional variable property laminar flow between two parallel plates. For 
this model, two-dimensional Navier-Stokes equations in fully conservative form are used to 
describe the flow field. These equations, in physical domain, can be written as [27,29] 


<9U dF dG ' 

aT + a7 + a7 + H = 0 


(i) 


where vectors U, F, G and H are written as 


P 

pu 

U = pv 

P E 

.pfj 

pu 

pu 2 + p + 

T XX 

F = puv -I- T xy 

(pE + p)u + t xx u + T zy v + q cx 4- <// 1 X 
pufj - pi 

pv 

PUV + Ty X 

G = PU 2 + P + Tyy 

(pE + p)v + T yy v Tx v u ^ Q c y 9/ty 
pvfj-pDty 
0 

0 

H = 0 

0 

-w, 


V- 
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The viscous stress terms appearing in the definitions of F and G are given in [27]. The relations 
for conduction heat transfer in x and y directions are given by 


, dT 

qc > = “ k a7 ' 


tiCy 


= -k 


dT 

dy 


(2) 


The terms qR X and qR y represent radiative fluxes in x and y directions, respectively; expressions 
for these are provided in the next section. The total energy flux in a given direction is given by 
the corresponding last term in the definitions of F or G. Consequently, this formulation involves 


all kinds of energy interaction including frictional (aerodynamic) heating. The coefficient 
of viscosity is evaluated by using the Sutherland’s formula and the coefficient of thermal 
conductivity is calculated by using a constant value of the Prandtl number. The total internal 
energy E appearing in U, F, and G is given by 

9,2 m 

E = P/f + — + X/Mi (3) 

£ i=l 


k i 


Equation (1) can be used to obtain solutions for all kinds of compressible flows. However, 
boundary conditions and numerical procedures for different flows are quite different. For 
supersonic flows the inflow conditions are specified and outflow conditions are obtained by 
extrapolation. The boundary conditions used along the surfaces are u=o, v=o, dP /dy = o, and 
T=T W . 

The governing equations and boundary conditions for the supersonic flow through a channel 
with a compression-expansion ramp is essentially the same as for the parallel plate geometry. 
However, a strong shock is produced at the compression comer and the flow becomes highly 
reacting from the beginning of the X 2 ~coordinate. (Fig. lc). 

The basic governing equations for chemically reacting compressible flow through a circular 
tube can be written as [29,30] 


d\J dV 1 0 , ,,, 

it + t + ~ y (; ) = H/y 

■ d t ax y ay 


( 4 ) 
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where x and y represent the stream wise and radial coordinates, respectively. In Fig. lb, the 
radial coordinate is denoted by r. Thus, both r and y notations are used to represent the radial 
coordinate for the circular geometry. The definitions of vectors U, F, and G in Eq. (4) are same 
as given in Eq. (1) and vector H is expressed as 

H = [0 0 p + rgg 0 y wj] (5) 


For the circular tube geometry, the viscous terms appearing in Eq. (1) are given by 


-->(! 


du u\ du 
+ dy + y) h dx 


. ( du du v\ dv 
~ \ dx + dy + y) ^ dy 


( du dv\ 

\dx dy yj y 


The boundary conditions for the circular tube geometry are similar to the parallel plate geometry. 


Radiative Transfer Models 


Evaluation of the energy equation presented in Eqs. (1) and (4) requires an appropriate 
expression for the net radiative flux in each direction. A suitable radiative transport model is 
needed to represent the true nature of participating species and transfer processes. In this section, 
a brief discussion of various absorption models is given and essential equations for the radiative 


flux are presented. 
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Absorption Models 

Several models are available in the literature to represent the absorption-emission character- 
istics of molecular species. The total band absorptance of a vibration-rotation band is given by 

fOO 

A = / [1 - exp£« w X)d(w - u; 0 )] (7) 

where k u is the volumetric absorption coefficient, u> is the wave number, u 0 is the wave number 
at the band center, X = Py is the pressure path length, and the limits of integration are over the 
entire band pass. Various models are used to obtain the relation for A in Eq. (7). 

The gray gas model is probably the simplest model to employ in radiative transfer analyses. 
In this model, the absorption coefficients is assumed to be independent of frequency, i. e., Km 
is not a function of w. A convenient model to represent the average absorption coefficient of a 
gray gas is the Planck mean absorption coefficient k p which is defined as [1] 

too 

Kp= K eUTI^MT) < 8a > 

J o 


For a multiband gaseous system, this is expressed as 

Kp = [P j /(aT , )]^ £ „.(T)S i (T) (8b) 

i 

where Pj is the partial pressure of jth species in a gaseous mixture, e^iCO is the Planck function 
evaluated at the ith band center, and Sj(T) is the integrated band intensity of the ith band. 

As defined in Eq. (8), « p is a property of the medium. When /c p is evaluated at the 
temperature of the gas, it is actually a mean emission coefficient and it becomes equal to the 
actual mean absorption coefficient only for the conditions of equilibrium radiation field. For a 
nonuniform temperature field, the mean absorption coefficient used for the optically thin radiation 
is the modified Planck mean absorption coefficient which for black bounding surfaces is defined 
as [1,9] 



K ul {T)e bul (T w )(L> 


Km(TT w ) = 


/db{T w ) 


(9a) 
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Note that K m is a function of both the gas temperature and the wall temperature. An approximate 
relation between /c p and is available for infrared radiation as [1] 

*m(T,T w ) = k p (T w )(T w /T) (9b) 

This expression is usually employed in gray gas radiative energy transfer analyses. 

Several other models for the mean absorption coefficient are available in the literature [1,31]. 
Since these models account for detailed spectral information of molecular bands, this approach of 
radiative formulation is referred to as the “pseudo-gray formulation.” The gray gas formulation 
for radiative transport is very useful in parametric studies. 

There are several nongray models available in the literature to represent the absorption- 
emission characteristics of vibration-rotation bands. These are classified generally in four classes, 
(1) line-by-line (LBL) models, (2) narrow band models, (3) wide band models, and (4) band 
model correlations. A complete discussion on usefulness and application of these models is 
provided in [12, 13], For many engineering applications, wide band model correlations provide 
quite accurate results. The most commonly used wide band model correlations are due to 
Edwards [5, 11] and Tien and Lowder [9]. The Tien and Lowder correlation for the total band 
absorptance is a continuous correlation and is given by the relation 

A(u J)= A(u,/?)/A 0 = ln{uf(t)[(u + 2)/(u + 2f(t))] + 1} (10) 

where 

f(t) = 2.94[1 - exp(— 2.60t)], t = £/2 

and u = SX/A 0 is the nondimensional path length, (3 = 2ny/d is the line structure parameter, 7 
is the line half width, S is the integrated band intensity, and A 0 is the band width parameter. 
Equation (10) provides accurate results for pressures higher than 0.5 atmosphere [12, 13]. 

Spectral properties and correlation quantities for various radiation participating species are 
available in [5, 9, 1 1]. These are useful in gray as well as nongray radiative formulations. 
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Radiative Flux Equations 

For many engineering and astrophysical applications, the radiative transfer equations are 
formulated assuming one-dimensional planar systems. For diffuse nonreflecting boundaries and 
in absence of scattering, the expression for the spectral radiative tiux in the normal direction 
is given by [1, 19] 


qRw(y) = eiu, - e 2 

+ 


-II 

- J V 


f y 

' 3 , 

/ FludZ^wt-Xp 

- y) 

*o 


f l , ^ 

3 , 

/ F^z)*^ exp 

-y) 

h 



(Iz 

dz 


( 11 ) 


where 


— ®oi (z) F'2 w(^) ^ 2 (z) *-2u; 


It should be pointed out that the exponential kernel approximation has been used in obtaining 
Eq. (1 1). The total radiative flux in a given direction is expressed as 


r 00 

<IR = / <tllu, du 


( 12 ) 


A combination of Eqs. (11) and (12) provides a proper form of total radiative flux equation 
for obtaining nongray solutions of molecular species. Any convenient absorption model can be 
used to obtain nongray results. 

For a gray medium, Eq. (11) reduces to a simpler form and upon differentiating the resulting 
equation twice, the integrals are eliminated and there is obtained a nonhomogeneous ordinary 
differential equation as [1, 16, 27 1 


1 d 2 q R (y) !) x :) de(y ) 

-7 — r-r- ~ 7 ( lit(.y) = j — 

K“ dy- I k dy 


(13) 


where k=k p . Equation (13) is a second order differential equation and, therefore, requires two 
boundary conditions. For nonblack diffuse surfaces, these are given as 



)[<iR(y)] v= o 


1 


d, 


lit 


dy 


v =0 


= 0 


(14a) 



10 


(H)(«WW -£[£]-» 


(14b) 


Equation (13) along with boundary conditions can be used to obtain the energy equation for 
gray gas radiative interaction. For black walls and Ti = T 2 , the boundary conditions for Eq. 


(13) become 


q R (l/2) = 0,^q R (0) = — (dq R /d £) f=0 > T ° = = y/L 

• To 


For a black circular tube, the spectral radiative heat flux in the radial direction is given by 


the expression [17] 


q ^ (r)= ?r{£./“ MK "“ xp [^ (r - 

— / F u ,(r / ) K^aexp — (r' — r) dr 

Ji cos 7 


; — r ( ) dr 


+ / F w (r') Ku»aexp — (r + r' — 2r sin 7 )dr| > d 7 

•/ rsin 7 L cos 7 J J 

where F u ,(r , )=e u ;(r / ) - eu;(Tw) and constants a and b have values of unity and 5/4, respectively. 
A combination of Eqs. (12) and (16) provides a convenient form of the total radiative flux for 
nongray analyses. 

For a gray medium, the expression for the total radiative flux can be obtained from differential 


approximation as [1, 17] 


arid, j 9, dx 4 

j; 7d? (r,|ll) _ t k " <w = 


For a black tube, the boundary conditions for Eq. (16) are found to be 

?q R (l) = -T- T-TT^qit) ' <llt(0) = 0 , To = n\'o,( = r/r 0 (18) 

2 r o Ls J £=i 

Equation (17) is used along with Eq. (18) for general one-dimensional gray gas formulation 
and analyses. 

With certain modifications, the radiative flux equations presented in this section can also 
be used to investigate the radiative interactions in the flow direction. The procedure for doing 
this is provided in [27], 
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Chemistry and Thermodynamic Models 


Chemical reaction rate expressions are usually determined by summing the contributions 
from each relevant reaction route (or path) to obtain the total rate of change of each species. 
Each path is governed by a law of mass action expression in which the rate constants can be 
determined from a temperature dependent Arrehenius expression. The reaction mechanism is 
expressed in a general form as 




j=l 


kbi 


E^Cj, 


j=i 


= 1, nr 


( 19 ) 


where ns= number of species and nr= number of reactions. The chemistry source terms in Eqs. 
(1) and (4) are obtained, on a mass basis, by multiplying the molar changes and corresponding 
molecular weight as 


nr 


w, = M,Cj = Mj ( 7 ," - 7 ij) 


1 = 1 


no 

k,. n c*r - hi n <#■ 


(20) 


m=l 


m=l 


The reaction constants kf t and k bi appearing in Eqs. (19) and (20) are determined from an 
Arrhenius rate equation given by 

h = Ai T Nl exp I 


<-§) 


(21) 


where 


^bi — kfj/keqi , k e qj 


-fe) 


AN 


exp 


•AGr; 

RT 


The coefficients A, N, and E appearing in Eq. (21) are given in Table 1 and the Gibbs energy 
term AGr, is calculated as 

ns ns 

AGr, = O'ijSi - 53 T, j& - J = 1 - nr (22) 

j=l j=l 


| = Aj(T - In T) + 5i T 2 + ^ T 3 


R 


Di 


+ST 4 + |T 5 + F j + GjT 
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20 


where 
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The gas constant for the mixture is evaluated by a mass-weighted summation over all species 
as 

R = E f ' R i < 23 > 

j=i 

The equation of state for the mixture is written as 


P = 4t 


(24) 


Method of Solution 


The governing equations are transformed from the physical domain (x, y) to a computational 
domain (£, 7?) using an algebraic grid generation technique similar to the one used by Smith and 
Weigel [32]. The grid spacing is kept uniform in the flow direction and compressed near the 
boundaries in the normal direction. The governing equations, Eqs. (1) and (4), are expressed 
respectively in the computational domain as 


d\) dF dG - n 
aT + di + ^ + H = 0 


(25) 


where 



dF 

dt 


ld_ 

^ r] dr] 


(r/G) = ft/. 


U = UJ, F = Fy n - G x v , G = Gx ( - Fy^, 
H = HJ , J = y,, - x. 


(26) 


The temporal descretization procedure for Eqs. (25) and (26) in given in [25-27]. Once 


this has been performed, the resulting system is spatially differenced using the explicit unsplit 


MacCormack predictor-corrector scheme [33]. This results in a spatially and temporally discrete, 
simultaneous system of equations at each grid point [34,35]. Each simultaneous system is solved 
using the Householder technique [36,37], and is marched in time until convergence is achieved. 
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The details of gray as well as nongray radiative flux formulations and solution procedure 
are available in [19,27]. For the gray gas model, the governing equations are discretized using 
a central difference scheme. The descretization of Eqs. (13) and (17) results, respectively, in 


Aj, -1 


2 

Bj 




1 \ ') , 
A yj + T fcj 




(28) 
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a] + -yjC, 


qj+i - 
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qj-1 = BBS 


(29) 


where 


Aj = Ay;(l + ft) , ft = A yj (l + ft) , Cj = ftAyj , 


RIIS = 1.5/cj 


e j+i 


- e; 


L ^jAyj 


+ 


e j c j-i 
Ayj 


, A yj = yj — yj— i, ft = 


w - yj 
yj - yj-i 


Equations (27) and (28) along with boundary conditions given by Eqs. (15) and (18) form 
tridiagonal systems of equations which can be solved efficiently by the Thomas algorithm. 


In the nongray gas formulation, the divergence of the radiative flux is evaluated using a 
central differencing scheme and is treated as radiative source term in the energy equation. Since 
the radiative flux terms involves integral formulation, unlike other flux terms which are only it] 


a differential form, it is uncoupled and treated separately. 


Physical Co nditions and Data Source 

The physical conditions for which specific flowfield analyses and computations are needed 
are discussed in [21-28]. In this work selected parametric studies have been conducted for 
certain flow and physical conditions. Radiation participating species considered are H 2 0, OH, 
and NO. Radiative properties of these species are available in [5, 9, 11-13]. Different amounts 
of these gases, in combination with air, are considered for parametric studies. Essential data for 
the chemistry model employed are obtained from Refs. 38-40 and these are provided in Table 1. 



14 


For basic studies, the physical dimensions considered for the channel are L=3 cm and L,=10 
cm, and for the circular tube they are L=D=3 cm and L»= 10 cm. Certain results, however, 
were obtained also for other dimensions. 

Results and Discussion 

Based on the theory and computational procedure described in the previous sections, an 
existing computer code was modified to solve the two-dimensional Navier-Stokes equations 
for radiating supersonic laminar flows between two parallel black plates. A similar code was 
developed for radiating supersonic flows in a circular tube. In most cases, the radiative interaction 
was considered only in the normal direction. Extensive results have been obtained for pure H 2 O, 
OH, and NO as homogeneous participating species, and for different mixtures of these species 
with air. Selected results are presented and discussed in this section. . 

For the parallel plate geometry (3 cm x 10 cm), a comparison of the divergence of radiative 
flux for general (nongray), gray, and their optically thin limit models is presented in Fig. 2 for 
two different y- locations (y=0.2 and 1.5 cm). The inflow conditions for this case are Poo = 1 atm, 

= 1,700 K, Moo = 4.3, f Hj o = 0.5, f 0j =0.1, and f Nj = 0.4. The gray formulation is based 
on the modified Planck mean absorption coefficient which accounts for the detailed information 
on different molecular bands. The range of optical thickness calculated in [27] was found to 
be between 0.0003 and 0.4. Thus, for the physical model and inflow conditions considered, 
the radiative interaction is essentially in the optically thin range. No significant difference in 
results is observed for the two y-locations. This is a typical characteristic of the optically thin 
radiation [1,10]. The solution of the gray formulation requires about ten times less computational 
resources in comparison to the solution of the nongray formulation [27]. 

As mentioned earlier, the Planck mean absorption coefficient /c p (or /c m ) is considered to 
be an optically thin radiation absorption coefficient, although it has been used in other optical 
ranges as well [1,9]. The appropriate absorption coefficient for the optically thick radiation is the 
Rosseland mean absorption coefficient kr. It has been pointed out in [1] that if the medium is 
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gray then kr=acp=k; otherwise /c p >/cr. Thus, use of « p (or K m ) in pseudo- gray gas formulation 
will provide maximum influence of radiative interaction. The k p values for H 2 O, OH, and NO 
have been calculated from Eq. (8b) by employing radiative properties available in [5, 9, 11] 
and these are illustrated in Fig. (3). Figure 3a shows the results of 100% homogeneous species 
whereas results in Fig. 3b are for different mixtures. The results provide indication of radiative 
ability of different species at a given temperature. Values of k p for other species are available 
in [1, 9], 

The results of supersonic entrance region flow between parallel plates are presented in Figs. 
4-10 for different physical and inflow conditions. As mentioned earlier, the basic physical 
dimensions considered for the channel are L=3 cm and L x =10 cm. In most cases, the inflow 
conditions considered are Poo=l atm, Too=l>700 K, Uoo=2574 m/s (Moo=3.0) with varying 
amounts of radiation participating species in combination with air. Certain variations in physical 
and inflow conditions are also considered for parametric studies. The chemical reactions are 
not considered in obtaining the results of Figs. 4-10, and radiative flux results are presented 
only for the normal direction. 

Results for radiative flux are illustrated in Figs. 4 and 5 as a function of the nondimensional 
y-coordinate. For P=1 atm, the results presented in Fig. 4 for different water vapor concentrations 
indicate that the radiative interaction increases slowly with an increase in the amount of the gas. 
The results for 50% H 2 0 arc presented in Fig. 5 for two different pressures (Poo=l and 3 atm) 
and x-locations (x=5 and 10 cm). It is noted that the increase in pressure has dramatic effects 
on the radiative interaction. 

For a mixture of 50% H 2 0 in air, the conduction and radiation heat transfer results are 
compared in Fig. 6 for Poo = 3 atm and for two different x-locations (x=5 and 10 cm). The results 
demonstrate that the conduction heat transfer is restricted to the region near the boundaries and 
does not change significantly from one x-location to another. The radiative interaction, however, 
is seen to be important everywhere in the channel, and this can have a significant influence on 


the entire flowfield. 
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Comparative results for 100% homogeneous species of H 2 O, NO and OH are illustrated in 
Fig. 7 at the exit plane (x=10 cm). As would be expected, the radiative contribution of H 2 O (with 
five bands) is significantly higher than NO and OH. Only the fundamental bands of NO and OH 
are considered in this study, and it is noted that NO is a better radiating gas in comparison to OH. 

For a mixture of 25% H 2 O in air, radiative flux for two different plate spacings (L=3 and 
6 cm) are illustrated in Fig. 8 for two x-locations (x=5 and 10 cm). The rate of radiative 
transfer is a strong function of the amount of the participating species and the pressure path 
length, PL. consequently, the results for the larger plate spacing indicate significantly higher 
radiative interactions. 

The effect of increased Mach number on the radiative transfer is illustrated in Fig. 9 for 
pure H 2 O and for a mixture of 50% H 2 O in air. The results shown are for the exit plane (x=10 
cm). At higher Mach number, the boundary layer is relatively thinner and the temperature in the 
boundary layer is significantly higher. This, in turn, results in higher rate of radiative transfer. 

For a mixture of 50% H 2 O in air, comparative results for the parallel plate channel and the 
circular tube are presented in Fig. 10 for two x-locations (x=5 and 10 cm). The results for the 
circular tube in general, exhibit the same trend as for the parallel plate geometry. Since the 
circular geometry provides additional degrees of freedom for radiative interactions [17], the rate 
of radiative transfer is higher for the tube. 

The influence of radiative interactions in chemically reacting supersonic internal flows 
was investigated by considering the physical model shown in Fig. lc. The specific problem 
considered is the supersonic flow of premixed hydrogen and air in a channel with a compression 
corner on the lower boundary. The physical dimensions considered for obtaining results are 
L=2 cm, xi= 1 cm, xj-l cm, Lx=xj+X 2=3 cm, and a=10°. The inlet conditions, which are 
representative of the scramjet operating condition, are P=1 atm, T=900 K and M=4.0. The 
flow is ignited by the shock from the compression comer. The flowfield for this problem has 
been investigated by several researchers [21-28] where different chemistry models have been 


used. Influence of radiative interactions was investigated in [25-27] by considering a simple two- 
step chemistry model. Recently [28], a comparative study of the flowfield was conducted by 
employing three chemistry models (see Table 1). It was found that significant amount of radiation 
participating species are produced by the 35-step chemistry model. The results presented in Figs. 
3-10 provided essential information on radiative behavior of these species. Complete discussions 
on use of the three chemistry models with and without radiation are provided in [28]. Selected 
results are presented here to demonstrate the influence of radiative interactions. 

The computed results for the 35-step chemistry model are presented in Figs. 1 1-14 with and 
without radiative interactions. The results were obtained by using a 31x31 grid; this was found 
to be an appropriate grid for the model. The variations in temperature, pressure, and species 
concentrations along the x-coordinate are shown for a y-location of 0.02 cm from the lower wall. 
The shock occurs at about x/L x =l/3, and it is noted that the 35-step chemistry model predicts 
the ignition time accurately. The temperature is seen to increase uniformly along the channel 
(Fig. 11) and there is a significant increase in pressure after the shock (Fig. 12). Results of 
Figs. 13 and 14 show that significant amounts of radiation participating species H 2 0 and OH are 
produced after the shock. The effect of radiative interaction is to lower the amount of species 
production due to radiative transfer in the x-direction. 

The results of radiative transfer by the three chemistry models are compared in Figs. 15 and 
16 at y=0.02 cm. The results for the normal radiative flux presented in Fig. 15 demonstrate that 
radiative interactions increase rapidly after the shock. The three models are seen to predict the 
same general trend. The results of streamwise radiative flux illustrated in Fig. 16 show that the 
net q R * decreases towards the end of the channel. This is due to cancellation of fluxes in the 
positive and negative x-directions. It is noted that the net radiative transfer is in the negative 
x-direction. The 18-step and 35-step models are seen to predict significantly higher qRx than 
the 2-step model. This is because radiative heat transfer is a strong function of temperature, 
pressure and species concentration which are higher (in the positive x-direction) for the 18-step 
and 35-step models than the 2-step model. 
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Conclusions 

Two-dimensional compressible Navier-Stokes equations have been used to investigate the 
influence of radiative energy transfer on the entrance region flow under supersonic flow condi- 
tions. Computational procedures have been developed to incorporate gray as well as nongray 
formulations for radiative flux in the general governing equations. Specific results have been 
obtained for different amounts of H 2 O, OH and NO in combination with air. Results demon- 
strate that the radiative interaction increases with an increase in pressure, temperature and the 
amount of participating species. This can have a significant influence on the overall energy 
transfer in the system. Most energy, however, is transferred by convection in the flow direction. 
The radiative interactions in reacting flows have been investigated by considering the supersonic 
flow of premixed hydrogen and air in a channel with a compression comer at the lower bound- 
ary. Depending upon the chemistry model employed, the radiative interaction is seen to change 
significantly the temperature, pressure and species concentration in the flow direction. 
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Table l. Hydrogen- Air Com bullion Mechanism [40) 


REACTION 

A(mole«) 

N(tm’) 

E(caloncs/gm-niGlc) 

following reactions commute the 18-step model M 


(1) Oj + ih^OH + OH 

1.70 x 10" 

0 

48150 

(2) Oj « II h OH + O 

1.42 x 10" 

0 

16400 

(3) H a ♦ Oil +■• H 2 0 ♦ H 

3.16 x I0 7 

18 

3030 

(4) Hi + 0 «-# OH ♦ H 

2.07 x 10" 

0 

13750 

(5) OH + OH 4- HiO ♦ O 

5.50 x I0' 1 

0 

7000 

(6) 11 + OH ♦ M ^ H a O ♦ M 

2 21 x I0 71 

—20 

0 

(7) Ilf li f M ♦ > Hi ♦ M 

6 53 x I0 ,J 

— 10 

0 

(8) Mo HOj t M 

3.20 a 10 11 

10 

0 

(9) Oil ► IIOs *+ Oj f ll 2 0 

S .OO x 10" 

0 

1000 

(10) H + HO i Hi + 0 2 

2.53 x 10" 

0 

700 

(11) H + HOj w OH ♦ OH 

1.99 x 10" 

0 

1800 

(12) O *■ HOj m Oj ♦ OH 

5 00 a 10 15 

0 

HXX) 

(13) IIOjtHOj t-*Oj + HjOj 

1.99 x 10" 

0 

0 

(14) H 2 -r HOi ^ H + HjOj 

3.01 x 10" 

0 

18700 

(15) OH + HiOa +♦ HaO ♦ HOj 

1.02 x 10" 

0 

1900 

(16) H + HiOj <-» H a O ♦ OH 

5.00 x 10" 

0 

10000 

(17) 0 + lliOa^>OH + HOi 

1.99 x 10" 

0 

5900 

(18) H 2 Oa + Mh OH + OH ♦ M 

1 21 x I0' 7 

0 

455(X) 

♦♦ remaining reactions complete the 35 step model •• 


(19) 0 2 + M«-*0 + 0 + M 

2.75 x 10” 

— 10 

118700 

(20) N 2 fMnNfN>M 

3.70 x 10“ 

— 16 

225000 

(21) N + Oa ^ O ♦ NO 

6.40 x 10* 

1.0 

63(X) 

(22) N ♦ NO O ♦ Na 

1 60 a 10 11 

0 

0 

(23) N + Oil hIU NO 

6.30 x 10" 

05 

0 

(24) 1U NO > M m UNO + M 

5 40 x 10" 

0 

— 6<X) 

(25) H > UNO «-» Hj ♦ NO 

4 80 x 10" 

0 

0 

(26) O ♦ UNO ~ OH ♦ NO 

5 00 x 10" 

05 

0 

(27) OH ♦ UNO ** HaO + NO 

3 60 x 10" 

0 

0 

(28) HOi ♦ UNO 4-> H 2 Oj ♦ NO 

2 00 x 10" 

0 

0 

(29) HOi ♦ NO ^ OH t NOj 

3 43 x 10" 

0 

- 260 

(30) 11 + N<)i ^ OH + NO 

3 50 x 10" 

0 

1500 

(31) O + N(>i ^ Oj + NO 

1.00 x 10" 

0 

6(X) 

(32) NOa + Me»0 + NO+M 

1.16 x 10" 

0 

66(XX) 

(33) M + OH + NO ** HNOa ♦ M 

5 60 x 10" 

0 

— 1700 

(34) MtOiltNO]H UNO) + M 

3 00 x 10" 

0 

— 38410 

(35) OH + IINOj «■» IIjO + NOj 

1 60 x 10" 

• 

0 

0 

♦♦ following reactions constitute the global 2 step 

model [24, 27, 39 J •• 

(l") ll 2 ♦ 2 OH 

11.4 x I0 47 

— 10.0 

4865 

(2") 2 OH + lla ^ 2 HjO 

2 50 x 10" 

— 13 0 

42500 








Fig. 2 Divergence of radivc flu* along the channel for gray 
and nongray model*, M_ ■ 4.3. 
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3a Variatioo of the Planch-mean abeorption coefficient 
for HfO, NO, and OH. 



Pig. 3b 


Variation of the Planck-mean abeorption coefficient for 
diffluent mixture* of HjO, NO, end OH. 
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Fig. 8 Radiative flux ve. y for two different plate spacing* (L * 3 
and 6 cm) at x * 5 and 10 cm, M. * 3.0. 
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Fig. 10 Radiative flux vi. y for channel and tube (L » D ■ 3 cm) 
at x ■ 5 and 10 on, 80% HjO, M_ * 3.0. 
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Fig. 9 Radiative flux ve. y for 50% and 100% HjO, x • 10 cm, 
M« * 3.0 and 4.5. 


Fig. 11 Temperature variation with x for reacting, and reacting 
and radiating flowe. 
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Fig. 13 Pressure variation with x for reacting* and reacting and 
radiating Howe. 



Fig. 13 Variation of HjO maaa fraction with i for reacting, and 
reacting and radiatingltdws. 
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Fig. 14 Variation of OH mass fraction with i for reacting, and 
reacting and radiating flows. 
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